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Abstract continued; 


a given radar resolution, the expected classification accuracy is shown to be 
dependent upon both the general soil moisture condition and also the geographical 
distribution of land-use (field-size distribution and dispersion of categories) 
and topographic relief. An analysis of cropland, urban, pasture/rangeland, and 
woodland subregions within the test site indicates that multi-temporal detection 
of relative soil moisture change is least sensitive to classification error resulting 
from scene complexity and topographic effects. 

The LOO m by 100 m radar resolution is found to yield the most robust classification 
results, and it is concluded that further degradation of image resolution should be 
implemented in post-detection processing when and where coarse resolution analysis 
is warranted. 
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A SIMULATION STUDY OF SCENE CONFUSION PACTORS IN 
SENSING SOIL MOISTURE FROM ORBITAL RADAR 

M. C. Dobson, S. Moezzi, F. T. Ulaby, and E. Roth 
Remote Sensing Laboratory 
University of Kansas Center for Research, Inc. 

Lawrence, Kansas 66045-2969 

ABSTRACT 

Simulated C-band radar imagery for a 124-km by 108-kro 
test site in eastern Kansas is used to classify soil 
moisture. Simulated radar resolutions are 100 m by 100 m, 1 
km by 1 km, and 3 km by 3 km; all images are processed with 
greater than 23 independent samples. The simulated radar 
operates at 4.75 GHz with HH polarization and over 7° to 17° 
angles of incidence. 

Distributions of actual near-surface soil moisture are 
established daily for a 23-day accounting period using a 
water budget model dependent upon precipitation, potential 
evaporation, crop-canopy cover, crop development stage, 
surface slope, antecedent soil moisture, and soil hydrologic 
properties. Within the 23-day period, three orbital radar 
overpasses are simulated roughly corresponding to generally 
moist, wet, and dry soil moisture conditions. The radar 
simulations are performed by a target/sensor interaction 
model dependent upon a terrain model, land-use 
classification, and near-surface soil moisture distribution. 
Rayleigh fading, layover, and shadow are accounted for by the 
model. For each overpass date and each radar resolution, the 
received power and range position of a given pixel is used to 
classify near-surface soil moisture via a generalized 


algorithm requiring no ancillary data about scene 
character ist ics . 

The accuracy of soil-moisture classification is 
evaluated for each single-date radar observation and also for 
multi-date detection of relative soil moisture change. In 
general, the results for single-date moisture detection show 
that 70% to 90% of cropland can be correctly classified to 
within +/- 20% of the true percent of field capacity. For a 
given radar resolution, the expected classification accuracy 
is shown to be dependent upon both the general soil moisture 
condition and also the geographical distribution of land-use 
(field-si 2 e distribution and dispersion of categories) and 
topographic relief. An analysis of cropland, urban, 
pasture/rangeland, and woodland subregions within the test 
site indicates that multi-temporal detection of relative soil 
moisture change is least sensitive to classification error 
resulting from scene complexity and topographic effects. 

The 100 m by 100 m radar resolution is found to yield 
the most robust classification results, and it is concluded 
that further degradation of image resolution should be 
implemented in post-detection processing when and where 
coarse resolution analysis is warranted. 


1.0 INTRODUCTION 
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Simulation tachniquaa have been employed to etudy the 
relationehip between spatial resolution and the accuracy at 
which soil moisture can be estimated from orbital C-band 
radar imagery [1,2] . These studies were based upon the 
land-use and crop-canopy-cover distributions present within a 
relatively small agricultural test site (18 km x 19 km) 
adjacent to the Kansas River in eastern Kansas. Image 
simulation techniques were used to generate 
synthetic-aperture radar (8AR) images at a frequency of 4.75 
OHz with HH polarization and with angles of incidence between 
7* and 22* from nadir. S.*.P images were produced at three 
different spatial resolutions: 20 m by 20 m with 12 looks, 

93 m by 100 m with 23 looks, and 1 km by 1 km with 230 looks. 
In addition, simulated real-aperture radar (RAR) Imagery was 
produced with a spatial resolution of 2.6 km x 3.1 km with 
363 looks. Analysis of these images demonstrated that for 
relatively flat agricultural portions of the test site about 
90% of the 20-m by 20-m pixel elements can be correctly 
classified to within +/- 20 percent of field capacity using a 
generalized soil moisture algorithm. In general, moisture 
classification accuracy was found to be greatest for coarser 
resolution Imagery due to the Increased number of looks; 
however, the results also showed a distinct 
classification-accuracy dependence on the complexity of the 
"true" soil moisture distribution and also upon the spatial 
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distribution of land-use elements within the test site. 

As a consequence, the current study is designed to 

examine further the effects of the spatial distribution of 

land-use categories, the agricultural field-size 

distribution, the crop-canopy mix, and the variability of 

local topographic relief on the soil-moisture classification 

accuracy achievable by various orbital radar resolutions at 

4.75 GHz, HH polarization, and angles of incidence from 7° to 

17°. An area of 124 km by 108 km, including most of the 

Lawrence, Kansas USGS quadrangle (1:250,000), serves as the 

test site. The area includes large subregions dominated by 

urban features, mixed cropland, rangeland and pasture, or 

deciduous woodland. Simulated radar imagery of this test 

site at resolutions of 100 m by 100 m, 1 km by 1 km, and 3 km 

by 3 km are used to classify soil moisture, which is 

subsequently compared to the input "true" soil moisture. 

» 

Classification accuracies of each radar resolution are 
compared for the whole test site and also for each of four 
subregions related to different mixtures of land-use. Since 
the number of processed looks for all resolutions is large (N 
>.23), the relative classification accuracies of each 
resolution should be only minimally biased by fading 
statistics . 

The dynamic behavior of each 100 m by 100 m grid cell 
within the simulation test site is modeled over a 23 day time 
period with respect to near surface soil moisture, crop 
canopy cover, crop stage-of -growth, and soil surface 
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roughness. The input parameters to this model include static 
conditions such as topography and soil association and also 
dynamic components consisting of cropping practices and daily 
meteorological conditions. The cropping parameters are based 
upon a stochastic treatment of average crop calendar, field 
size distribution, and crop development while the 
meteorological data includes daily rainfall and potential 
evaporation. The output of this model consists of daily 
updates of near surface (0-5 cm) soil moisture and radar 
backscatter category which is approximately equivalent to a 
Level III land-use category [3). The model is run for a 23 
day period and the outputs are saved on 3 dates corresponding 
to hypothetical orbital overpasses each nine days apart. The 
overpass dates were selected independent of any consideration 
of 'jT!;'.tal mechanics but rather to represent three 
distinctive soil moisture distributions over the test site: 
very wet, moist, and dry. The above moisture classifications 
are very general, however, since the large size of the data 
base and the late spring time frame of the simulations leads 
to highly variable regional soil moisture distributions on 
any given date. 

For each orbital overpass, a target-sensor interaction 
model produces simulated radar imagery for each of the three 
radar resolutions. The simulation model accounts for the 
effects of Rayleigh fading and geometric properties such as 
layover and shadowing [22]. Each simulated radar image is 
then subjected to a generalized algorithm (requiring only the 


amplitude of received power and the range position of a given 
pixel) which classifies the image into estimated soil 
moisture. These distributions of estimated soil moisture are 
subsequently compared with the distributions of actual 
near-surface soil moisture on a grid-cell basis for each 
date. 

In addition to testing the absolute classification 
accuracies of each radar resolution for each of the three 
overpass dates in an instantaneous sense, multi -temporal data 
from two of the overpasses is used to evaluate the merits of 
relative change detection of near surface soil moisture as 
estimated from each of the three simulated radar resolutions. 
The above process is shown schematically in Figure 1. 


2 . 0 TEST-SITE DATA BASE 


In order to quantify the radar backscatter ing from a 
given terrain element, certain geometric and dielectric 
properties of the target scene must be known. First, the 
three-dimensional cartographic coordinates of each element 
must be specified relative to the orbital radar in order to 
compute range, area, and local incidence angle. Secondly, 
the radar backscatter ing category must be established; this 
is roughly equivalent to a level-III land-use classification 
category [3J. Finally, many land-use categories have 
backscatter ing propertiis that vary as a function of 
crop-canopy cover, row directionality, and near-surface soil 
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Figure 1. Flowchart of simulation approach 












moisture, A three-tiered digital data base is constructed to 
describe the spatial distribution of category elements and a 
dynamic model acts upon this distribution to vary target 
dielectric and backscatter ing properties as a function of 
time. It is assumed that all target properties are laterally 
homogeneous within a given 100 m by 100 m terrain element. 

2 . 1 Terrain Model and Radar Backscatter ing Categories 

Digital elevation data from the Defense Mapping Agency 
provide a static model of the terrain geometry. These data 
are corrected for scanning errors and resampled to yield a 
mean elevation for each 100-m by 100-m grid element within 
the 124-km by 108-km test 3ite. An image-format presentation 
of the digital elevation data is shown in Figure 2. 

The specification of radar backscatter ing category for 
each 100-m by 100-m grid element involves a three-step 
process that accurately describes the spatial distribution of 
the categories shown in Table 1 in a stochastic sense. A 
two-dimensional digital matrix of Level-II land-use 
classification is given by USGS land-use and land-cover 
digital data (LUDA) for the Lawrence, Kansas quadrangle. 
Level-II categories with similar radar backscatter ing 
properties (such as lakes and rivers) are redefined as 
equivalent backscatter ing categories. The Level-II LCJDA 
category of cropland is insufficient to specify unique 
backscatter ing characteristics; thus a stochastic process is 




TABLE I. Area Percent of Data Base Assigned to 
Each Target Class 





































used to further define the spatial distribution of particular 
agricultural crops. A random sample of U-2 high-altitude 
color IR images is used to generate statistics on 
agricultural field-size distribution for each of the twelve 
counties within the test site. These statistics are then 
used to assign random f ie Id -boundary networks within each 
county. The distribution of field sizes is given by county 
in Table 2. 

Specific crop categories and row directions are randomly 
assigned to each field within a county, based upon an 
historical enumeration of crop acreage for each county 
provided by the Kansas State Board of Agriculture and the 
Missouri Department of Agriculture. These acreages are given 
by county in Table 3. In addition, since all crops are not 
grown concurrently, crop calendar data [4J is used to factor 
planting and harvest into the time history of each field. 
Within a given crop, planting and crop-development stages 
established for this area are used to change a given field's 
backscatter category from bare soil to that of the crop after 
emergence in a stochasitc fashion. The fields of each crop 
type are subdivided into ten subgroups each with a 
distinctive cropping history. Thus, the crop-type 
distribution will vary locally as a function of time within 
the 23 day simulation period. The land-use and crop-type 
distributions for the entire 124 km by 108 km test site are 
shown in Table 1 for each of the hypothetical orbital 
overpass dates. The simulation period runs from May 18 
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TABLE 2. Field-Size Distributions for the Agricultural Portions 
of the Land-Use Subregions 


Percent of Agricultural Area 


Field Size in Acres 


Subregion 

10 

20 

30 

40 

60 

80 

100 



120 

140 

160 

Urban (Kansas City) 

20 

18 

10 

15 

n 

16 

3 

2 

2 

7 

Pasture/Rangeland 

D 

11 

6 

18 

8 

- — i 

23 

3 

-J 

5 ! 

3 

14 

Cropland 

20 

23 i 

12 

19 

a 

6 

2 ; 

2 ; 

2 i 

3 

Woodland 

20 

23 

12 
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TABLE 3. Relative Percent of County Cropland 
Devoted to a Given Crop or Pasture/Range [9,10] 


Group A ® Anderson County 

Group B = Bates, Douglas, Franklin, Linn, and Miami Counties 
Gropu C = Cass, Jackson, and Johnson Counties 


Percent of Total Agricultural Land 


Group 

■ Wheat 

Sorghum 

Com 

Oats 

Soybeans 

Alfalfa 

Pasture 
Hay & Range 

A 

8.6 

7.4 

5.7 

0.5 

21.1 

13.7 

43.0 

B 

6.3 

9.6 

5.5 

0.5 

15.2 

10.3 

52.0 

C 

4.1 

5.2 

5.8 

1 

0.4 

11.6 

8.2 

64.7 


Note: Urban Subregion consists of most of Jackson and Johnson Counties 

Cropland Subregion consists of parts of Douglas, Franklin, Johnson, 
and Miami Counties 

Pasture/Rangeland Subregion consists of most, of Cass County 
Woodland Subregion consists of a large part of Linn County 
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(Julian day 138) until June 9 (Julian day 160) during which 
time corn and soybeans are emerging and this is reflected in 
Table 1. Examples of land-use and crop-category distribution 
are shown in Figure 3 for Julian day 141. 

2 . 2 Dynamics of Soil Moisture Distribution 

The above two components of the data base define the 
geometric properties of the test site and the distribution of 
backscatter ing categories. In addition, it is necessary to 

model certain dynamic conditions that largely determine the 
dielectric properties of the scene elements. Of major 
importance is the near-surface soil moisture of each 100-m by 
100-m pixel element as a function of time. 

The soil moisture is governed by soil type, local slope, 
crop canopy cover and stage of growth, antecedent soil 
moisture, precipitation, and potential evaporation. The 
distribution of soil types as generalized by soil 
associations from USDA/SCS county soil surveys is shown in 
Figure 4. The local crop calendar is derived for this area 
from historical records [4] and used to establish the daily 
transpiration rate for a given crop. Daily weather records 
from each of 25 reporting stations are used to generate 
digital overlays of daily precipitation (Figure 5) and 
potential evaporation. A water-budget model is used to 
update near-surface soil moisture on a daily basis for each 
grid cell. Finally, a normally distributed random-noise 
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Land-use: urban features are bright while water and woods 

are darkest. 


Figure 3. Land-use and crop-category distributions on Julian day 141 
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Enlargement of upper-left corner shows 51.2 km by 5! 
of total scene. 


Enlargement of 2b shows Kansas River and trees as 
black, urban features as white, the remainder of 
the image shows cropland of which soybeans are 
emphasized to show the presence of both north-south 
and east-west row directions. 
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Figure 5. Image presentation of the areal distribution of 

rainfall within the test site on each Julian date 
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component is added to the modeled aoil moiature in order to 
eimulate local, within-field variance in true aoil moiature 
£ 2 ]. The details of the aoil water accounting model and a 
listing of the computer program are given in Appendix A. 

Examples of the 0-5 cm soil moiature distributions 
produced by the model are shown in Figure 6 for Julian days 
141, 150, and 160 in image format. The corresponding 
cumulative areal uiatr ibutiona are shown in Figure 7a for 
each date. The influence of crop cover on soil moisture 
distribution is shown in Figure 7b for Julian day 150. These 
distributions when combined with the terrain model and the 
spatial distribution of radar backscatter categories 
collectively drive the radar image simulations discu33ed 
below. 


3.0 RADAR IMAGE GENERATION 


The average return power reradiated from each 
laterally homogeneous grid ce-' is given by the radar 
equation 


P„ G 2 X 1 o° A 


<4it) 3 R" 


( 1 ) 


2 

where is the average transmitted power, G is the two-way 
antenna gain, X i3 the wavelength, o° is the radar cross 
section per unit area, A is the grid-cell area, and R is the 
range. For a given sensor configuration, P T , G, and X are 
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(a) Julian day 141 . 


Figure 6. 


Distribution of 
Black represents 


0-5 cm soil moisture across the 
undefined (zero) soil moisture. 


test sit 





(b) Julian day 150. 



(c) Julian day 160. 
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PERCENT OF FIELD CAPACITY 

(a) Net distribution for each date. 


Figure 7. Cumulative distributions of soil moisture on satellite 
overpass dates. 
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Day 150 All Crops 




constant. For each grid cell element, area and range are 
determined from the static terrain model. In addition, a° 
varies as a function of local angle of incidence, 
backscatter ing category, and near-surface soil moisture; and 
for the purpose of radar simulation, o° is given by empirical 
fits to experimental airborne and truck-mounted scatterometer 
data [6]. Examples of empirical radar backscatter dependence 
on target category, incidence angle and near-surface soil 
moisture are given in Table 4. Radar backscatter ing 
coefficient o° is shown graphically in Figure 8 as a function 
of local incidence angle 6 for selected categories and soil 
moisture conditions. 

The power actually received at the antenna P f is 
dependent upon signal fading and atmospheric scattering and 
adsorption. At 4.75 GHz the atmospheric losses are assumed 
to be negligible for most conditions. In addition, signal 
fading is assumed to be ^-square distributed with 2 N degrees 
of freedom where N is the number of independent samples for a 
given range and azimuth radar resolution [7]. Hence, 

p Y (2) 

r \2N / 

where Y is a random variable with x-squared distribution and 
2 N degrees of freedom. 

The radar image simulation model accounts for the 
geometric effects of layover and shadowing. Examples of 
simulated orbital radar imagery are shown in Figure 9 for the 
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Figure 8. Radar backscattering c at 4.75 GHz with HH polarization 

as a function of local incidence angle for selected moisture 
conditions. 
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soil moisture distribution present on Julian day 141 (Figure 
5a) at radar resolutions of 100 m by 100 m, 1 km by 1 km and 
3 km by 3 km. These images are ground-range presentations 
and is scaled in dB to facilitate the presentation of the 
large dynamic range in P r across the image swath (“ 48 dB). 
The radar illumination is from the west (left side of 
images). Due to the relatively steep incidence angles (7 W - 
17°), the angular decay in P^ is readily apparent across the 
swath from left to right. In general, areas of higher 
near-surface soil moisture as related to antecendent 
precipitation appear brighter on the images, and this is mo3t 
apparent as diagonal stripes related to storm tracks. Also, 
areas of tree canopy cover and water bodies tend to be dark 
on the imagery simulated for Julian day 141, while urban 
features tend to appear bright and are especially noticeable 
in the far range (right side of images). 

The simulated orbital imagery for the three radar 
resolutions are also shown in Figures 10 and 11 for Julian 
days 150 and 160, respectively. Julian day 150 represents 
the wettest overall soil moisture conditions as indicated in 
Figure 7, and hence the images appear brighter than those for 
Julian day 141 (Figure 9). In contrast, Julian day 160 is 
shown by Figure 7 to represent the driest overall soil 
moisture conditions, and thus the images in Figure 11 appear 
darker than those for Julian day 141 (Figure 9). 

It should be noted that for all of the above simulated 
images (Figures 9, 10, 11), the number of independent looks 
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(a) 100 m by 100 m radar resolution. 


Figure 10. Simulated radar imagery of the test site on Julian 
day 150. 


3 






(b) 1 km by 1 km radar resolution. 
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(a) 100 in by 100 m radar resolution. 


Figure 11. 


Simulated radar imagery of 


the test site on Julian day 160. 
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ia large (N 23). Hence, the variance in within a given 
portion of the scene is only minimally dependent upon signal 
fading and is mostly the result of variance in local 
topographic relief, radar backscatter category, and 
near-surface soil moisture. In a visual sense, the 
interaction of relief, category, and moisture yield quite 
different spatial patterns of P £ on each of the three 
simulation dates. This is best seen in the 100 m by 100 m 
radar resolution imagery. Figure 12 shows enlargements of 
the northwest (upper-left) quadrant of the 100 m by 100 m 
imagery for each of the three overpass dates. This quadrant 
encompasses the test site used in previous orbital radar 
simulations [1, 2, and 6]. These images illustrate the 
following: 

1) For nearly uniform soil moisture conditions, the 
variance in P f is dominated by local topographic relief and 
radar backscatter category. This condition is most closely 
approximated by Julian day 150 in Figure 12b. 

2) For variable soil moisture conditions, the scene 
variance in P r is most closely related to local soil moisture 
and radar backscatter category which tends to mask variance 
in P r related to local topographic effects. This condition 
is best seen on Julian day 160 (Figure 12c) since an extended 
period of evapotranspirative losses in soil moisture has 
enhanced the relative difference in P r from each radar 
backscatter category. 

The above indicates the potential for achieving certain 
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mapping objectives not rigorously addressed within the 
confines of this study. First, the potential exists to 
classify soil type within relatively flat agricultural 
portions of the test site from imagery acquired shortly after 
a nearly uniform and saturating rainfall event. In this 
case, near-surface soil moisture is high and largely 
controlled by soil hydraulic properties related to soil type. 
In addition, for high moisture conditions, the relative 
uncertainty in P f related to crop-canopy attenuation and 
canopy backscatter is expected to be small [5). Secondly, 
the potential for crop discrimination from orbital radar 
imagery can be expected to maximize (for this frequency and 
angle of incidence) when the differential evapotranspirative 
dry-down of each crop has enhanced the inter -crop variance in 
P r . This condition would exist five or more day9 after a 
rainfall event. 

4.0 SOIL MOISTURE CLASSIFICATION 

In order to classify soil moisture using the simulated 
radar imagery, a generalized soil-moisture algorithm is 
derived from all experimental data for bare and 
vegetation-covered soil conditions (excluding woodlands). 

The classification algorithm relates estimated soil moisture 
“fs to received power ? r as a function of incidence angle 9. 

M f s ° [ p r -a <0 )]/e(0) (3) 
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where 


a(6) = 9.67 + 0.840 - 4.59 x 1O _2 0 2 + 8.27 x 10"V, and 
3(0) = 0.161 + 9.38 x 10~S - 4.97 x 10 _ S 2 + 1.21 x 1O _5 0 3 . 

In this case, e is estimated from the range position of a 
pixel on the radar image, assuming spherical earth geometry 
and a constant mean elevation of the test site above sea 
level. Thus, the classification algorithm is "blind" with 
respect to true local incidence angle and to the actual 
backscatter ing category of any given pixel [6]. Application 
of this algorithm to the received power images yields maps of 
estimated soil moisture, an example of which is shown in 
Figure 13 for a radar resolution of 100 m by 100 m on Julian 
day 141. 

Given the above algorithm, orbital radar imagery can be 
used to classify soil moisture in two ways. First, the 
imagery obtained at any given radar resolutin on any single 
overpass date can be pssed through the general algorithm 
(Equation 3) to yield estimates of the absolute soil moisture 
distribution for that date. The second approach is to make 
use of the multi-temporal coverage provided by an orbital 
system to yield estimates of the relative change in soil 
moisture. The radiometric and geometric stability of the 
Seasat-A L-band imaging radar has shown that such a procedure 
is feasible and relatively uncomplicated from the standpoint 
of image registration [8]. The two approaches are not 
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Figure 13. Near-surface soil moisture as estimated for Julian day 141 
from simulated radar imagery with a 100 m by 100 m resoluti< 


mutually exclusive and both will be explored in the ensuing 
sections with respect to soil moisture classification error 
as a function of radar resolution and the geographic 
distributions of local relief and backscatter category. 

4 . 1 Single Date Soil-Moisture Classification Accuracy 

The accuracy of soil-moisture classification is examined 
by evaluating the difference between the true soil moisture 

A 

M~ and the estimated soil moisture M* . This is 
is f s 

accomplished through registration of the two images (such as 
Figures 6 and 13) and computation of the difference. Due to 
the geometric distortion inherent in the radar image-forming 
process, image registration by simple coordinate translation 
is only accurate to within about +/- 1.3 pixels (130 meters), 
and this registration error is proportional to changes in 
local elevation across the image swath. Hence, a procedural 
error is introduced into the comparative process which is not 
related to true classification error. Also, the magnitude of 
this procedural error is proportional to the local variance 
in the "true" soil moisture distributions as shown in Figure 
6 . 

In order to examine the effects of various land-use and 
field-size distributions, four subregions are identified 
within the test site and relate to an urban area, mixed 
cropland, pasture and rangeland, or woodland. Figure 2 shows 
the spatial locations of these subregions, and their land-use 



and field-size distributions are tabulated in Tables 1 and 2, 
respectively. All subregions contain more than 30% pasture, 
grass, and rangeland, and are distinctive primarily in terms 
of the percent area occupied by cultural features 
(residential, buildings, and roads), water, woodland, and 
crops. In addition, the rangeland/pasture subregion is 
characterized by a greater percentage of large fields as 
compared to the other subregions. Finally, Figure 2 shows 
that the woodland and the rangeland/pasture subregions are 
located in areas of relatively large local relief. 

An example of soil-moisture classification error is 
shown in Figure 14 for the 100-m by 100-m resolution radar on 
Julian day 141. Classification error E m is defined by 


where 

M^ g = true soil moisture, and 

A 

Mfs “ estimated soil moisture. 

Figure 14a shows the category classification map for the 
woodland subregion where wooded areas are black, water is 
dark gray, cultural features are white, and agricultural land 
and pasture/rangeland are generally light gray. The 
difference between actual soil moisture M fg and classified 

A 

soil moisture M^ g is mapped in Figure 14b. is linearly 

represented by graytone and thus, dark and white areas 
represent overestimation and underestimation of soil 
moisture, respectively. The large P from cultural features 




(a) Backscatter category map: woods are black, water bodies are dark gray 

cultural features are white, and agricultural areas are light gray. 


(b) Soil moisture estimate error E,,,: overestimates of soil moisture are 

dark, underestimates of soil moisture are white, areas with small 
estimate errors are gray. 


Figure 14. Soil moisture classification error ^ on Julian day 141 within the 
woodland subregion resulting from use of the "blind" classifier on 
100 m by 100 m radar imagery. 





at A leads to an over estimation of moisture, while the 
low ? r from woodland at B and water at C yields a low 
estimate of soil moisture. Median gray tones in Pigure 14b 
relate to small estimate errors. A comparison of Pigures 14a 
and 14b shows moisture-estimate errors to be highly 
correlated with the spatial location of specific land-use 
categories, especially cultural features, trees, and water. 
Image registration errors yield white or black rings around 
specific features. Hence, the spatial organization of such 
confusion categories largely determines the moisture 
classification accuracy of a given radar resolution for a 
given geographic land-use setting. 

The single date soil moisture classification error can 
be examined as a function of radar resolution, general soil 
moiture condition (overpass dat e) , and geographic subregion. 

The soil moisture classification error resulting from 
radar resolutions of 100 m by 100 m, 1 km by 1 km, and 3 km 
by 3 km is shown for the entire 124-km by 108-km test site on 
each overpass date in Pigure 15. For all general soil 
moisture conditions (overpass dates), the distributions of 
resulting from classification of the 100-m resolution 
imagery are more peaked and yet have longer tails than the 
corresponding distributions of E m for the coarser 
resolutions. These long tails are related to the presence of 
confusion categories such as urban features, woodland, and 
water. The effects of these confusion categories at the 
coarser resolutions (1 km by 1 km and 3 km by 3 km) are to 
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Figure 15. Soil moisture classification error resulting from 

each radar resolution for all moisture dependent pixels 
in the test site (excluding woods) on a) Julian day 141, 
b) Julian day 150, and c) Julian day 160. 





Percent Area 




broaden the error distribution. 

The tendency of the coarser resolutions (1 km by 1 km 

and 3 km by 3 km) to yield bimodal distributions of E in 

m 

Figure 15 with a secondary peak ranging from -30% to -50% is 

primarily related to the presence of cultural features which 

have large scattering cross-sections relative to agricultural 

and rangeland areas. These overestimates of local soil 

moisture result from averaging the large P from cultural 

targets over a larger area. Hence, the magnitude of this 

secondary peak is proportional to both the net area occupied 

by cultural features and the dispersion of such features 

within the total scene, and the size of E at this peak is 

m 

proportional to the ratio of cultural to P r agricultural. 

The associated absolute moisture classification 
accuracies of the three radar resolutions are shown in Figure 
16. In general, the 100-m by 100-m resolution is shown to 
yield the most accurate estimates of soil moisture. For 
example, use of the "blind" generalized moisture algorithm on 

A 

Julian day 141 yields M fg within +/- 20% of true 
moisture M fg for 68% of the area using a radar with a 100-m 
resolution, while only 60% and 58% of the area is classified 
within this error limit using radar resolutions of 1 km and 3 
km, respectively. In Figures 14 and 15, this result is shown 
to be related to the spatial distribution of land-use 
confusion categories. 

The differences in absolute classification accuracy 
between the three radar resolutions are also dependent upon 












general soil moisture condition. This effect is seen by 
comparing the results achieved for different overpass dates 
in Figure 16a, b, and c for overpasses on Julian days 141, 
150, and 160, respectively. The classification accuracy of 
the coarse resolution sensors (1 km by 1 km and 3 km by 3 km) 
is seen in Figure 16c to be significantly reduced relative to 
the classification accuracy achieved with the 100 m 
resolution radar. The local variance in true soil moisture 
M^ a and local received power P r are seen to be greatest on 
Julian day 160 in Figures 6 and 12, respectively. As 
previously stated, this is largely the result of the 
differential evapotranspirative dry-down rates of the various 
crop canopies constituting the scene. Thus, the within-scene 
variance in soil moisture M fg is highly correlated with the 
crop distribution given in Table 1 which is dispersed in 
agregates given by the field size distribution (Table 2). 
Hence, at radar resolutions coarser than field size a serious 
degradation in moisture classification accuracy can be 
expected for imagery acquired during periods of protracted 
evapotranspirative loss. 

The effects on moisture classification error of varying 
the local distribution of land-use confusion categories are 
demonstrated by comparing the error distributions for the 
four land-use subregions. The error distributions for the 
urban, pasture/rangeland, cropland, and woodland subregions 
are compared in Figure 17, based upon the 100-m resolution 
radar imagery for Julian day 141. When the error 


(b) Moisture Dependent Pixels Only 
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distribution is plotted for all 100-m by 100~m grid cells 
within each region (Figure 17a), large overestimates of 
moisture, primarily in the urban and woodland subregions, are 
related to the presence of cult al features such as 
buildings and roads and also related to the presence of water 
bodies since for these categories soil moisture is undefined 
and any moisture estimate for these categories is therefore 
an overestimate. In a similar fashion, large underestimates 
of M^ g , best exemplified by the woodland subregion, ate 
largely related to the presence of deciduous trees, which are 
assumed to fully attenuate backscatter ing from the soil at 
4.75 GHz. 

The exclusion of nonagr icultural categories (cultural 
features, water, and woodland) from the grid-cell comparisons 
of M fg to M fs yields highly peaked distributions centered 
around « 0 error as shown for each subregion in Figure 17b. 
The woodland still exhibits a larger area where soil moisture 
is underestimated than the other subregions and this is 
largely the result of locally saturated to flooded soil 
moisture conditions. The radar backscatter model treats 
fully saturated soil as a near specular surface similar to a 
water body, and hence P £ is low at off nadir indicence 
angles. As a consequence, soil moisture M fg is generally 
underestimated. Similar results are obtained for the other 
two overpass dates. 

The absolute classification accuracy for Julian day 141 
within each of the four land-use subregions is shown in 


Figures 18a and 18b from simulated radar resolutions of 100 m 
by 100 m and 1 km by 1 km respectively. As expected from the 
above and from the distributions of land-use categories and 
field-size given in Tables 1 and 2, Figure 18a shows that the 
greatest classification accuracy is achieved for the cropland 
subregion and the poorest for the woodland subregion. Based 
upon land-use and field-size dist:: ihut ior.s alone, one would 
expect a greater absolute classification accuracy for the 
pasture/rangeland subregion than for the urban subregion in 
Figure 18a; however, the greater local topographic variation 
present within the pasture/rangeland subregion (Figure 1) 
leads to moisture classification errors related to the 
variance in local slope, which is unknown to the "blind" 
classification algorithm. This same effect also suppresses 
the absolute classification accuracy for the woodland 
subregion which is also "hilly" in nature. 

For a 1-km by 1-km resolution radar, the combined 
effects of the spatial distribution of land-use categories 
(the relative mix of categories and their respective size 
distributions) and topographic relief upon absolute 
classification accuracy yield the results shown in Figure 
18b. For areas where local topographic relief varies over 
spatial dimensions of hundreds of meters, the 1-km by 1-km 
radar resolution will tend to average local slope-related 
variance in P r , and thus yield absolute classification 
accuracies greater than those achieved by a finer resolution 
sensor (such as 100 m by 100 m) . This appears to be the case 
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for the pastur e/rangeland and woodland subregions of the test 
site. For example, at an absolute accuracy level of +/- 20% 
of field capacity, Table 5 shows that the percent area 
correctly classified within this limit from the 100-m 
resolution radar is 71.3% and 64.1% for the pasture/rangeland 
and woodland subregions, respectively; and the percent area 
correctly classified from the 1-km resolution radar increases 
to 79.4% and 73.3% for the two subregions, respectively. 

Conversely, for areas character ized by a large number of 
dj. .persed cultural targets (with generally large P f ), the use 
of a coarse-resolution radar, such as 1 km by 1 km, is shown 
to degrade absolute moisture classification accuracy relative 
to that achievable by a 100-m by 100-m resolution sensor; 
this effect is demonstrated by the urban and cropland 
subregions. For example, in Figure 18 the effect of 
dispersed cultural features and field size distribution leads 
to a decrease in percent of the urban subregion which is 
correctly classified to within +/~ 20% of field capacity from 
77.9% (100-m radar resolution) to 70.1% (1-km radar 
resolution). In a similar fashion, the percent area 
correctly classified to within +/- 20% of field capacity for 
the cropland subregion decreases from 82.0% to 75.6% for the 
100-m and 1-km radar resolutions, respectively. 

The above results for Julian day 141 are not independent 
of general soil moisture condition and the spatial 
variability of soil moisture. The absolute soil moisture 
classification accuracies for each of the four subregions are 
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TABLE 5. Percent of Moisture Variant Area Correctly 
Classified to Within +/-20 of True Soil Moisture 
( | E | < 2QZ) 

' m 1 — 


Julian Day 


141 



150 



160 


Radar Resolution 

100 in 

1 km 

3 km 

100 m 

1 km 

3 km 

100 m 

1 km 

3 km 

Subregions 


■ 








Cropland 

£2,0 


74.3 

87.5 

91.0 

91.8 

77.4 

73.9 

75.6 

Urban 

77.9 

BO 

65.0 

75.3 

83.6 

86.9 

76.2 

63.7 

58.5 

Rangeland/Pasture 

71.3 

79.4 

80.4 

72.5 

80.4 

82.3 

77.0 

74.6 

77 A 

Woodland 

64.1 

73.3 

72.7 

72.6 

83.6 

86.0 

68.3 

68.5 

68.4 





















shown for Julian da/s 150 and 160 in Pigures 19 and 20, 
respectively . In addition, the percent area correctly 
classified to within +/- 20% of field capacity are also given 
for each dat*?> and radar resolution in Table 5. 

Por the generally wet soil conditions prevalent on 
Julian day 150, comparison of the results shown in Pigure 19 
and Table 5 as a function of orbital radar resolution 
indicates that estimate accuracy increases with the 
additional spatial averaging provided by the coarse 
resolution radars for all subregions. T^is is explained by 
the distribution of soil moisture for this date which is 
primarily governed by the antecedent rainfall pattern. Since 
a large quantity of rain fell within most of the test site 
just prior to the simulated orbital overpass, the local 
properties of slope, soil texture, and crop canopy condition 
have not had sufficient time to exert a large influence and 
vary local soil moisture distributions. As a result, the 
added spatial averaging provided by the coarser radar 
resolutions acts to increase classification accuracy by 
averaging small spatial scale noise effects related to local 
relief and variance in local radar backscatter ing category. 
This is true even for the urban scene; since at very high 
soil moisture conditions, the P r from wet agricultural fields 
approaches that from the cultural features. 

Within the four subregions, the dependence of soil 
moisture classification accuracy upon radar resolution is 
shown in Pigure 20 and Table 5 for the generally dry and 
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Spatially variable so; L moisture conditions prevalent on 
Julian day 160. Classification accuracy is shown to be 
independent of radar resolution for absoluv.^ i^'timate error 
lM^ g - less than 20% of field capacity for the cropland, 

rangeland/pasture, and woodland subregions. However, as 
radar resolution is degraded the areal percentage of the 
cropland and woodland subregions with large absolute estimate 


errors, IM fa - M fs l >_ 30, does increase significantly. This 


4 f 3 


■ is attributed to the large local variance in true soil 
moisture M fg within these subregions on Julian day 1.60. The 
most extreme example of local variance in M^ g is given by the 
urban subregion which exhibits a pronounced decrease in 
classification accuracy as radar resolution is degraded. 


4.2 Multidate Change Detection of Soil Moisture 

The preceeding section shows that absolute moisture 
classification accuracy from a single date orbital radar 
observation is limited by the presence of scene confusion 
factors within the imagery and their size and spatial 
dispersion relative to the radar resolution. Within the 
present discussion, scene confusion factors are defined both 
as the presence of scene elements for which soil moisture is 
unidentified such as buildings, roads, water bodies, etc. and 
also the occurance of variability in P r (0) from scene 
elements possessing equivalent soil moisture. The latter 
results from natural variability in topographic slope, crop 
canopy type and stage of growth, row direction, and surface 


roughness . 

In single date sensing and classification of soil 
moisture, the confusion effects of cultural features and 
water bodies can be minimized (but not eliminated) by spatial 
filtering. Two approaches are feasible. First a simple 
intensity slice of the received power P f (0) could be used to 
roughly define water (dark) and point targets such as 
buildings (bright) within the image, the remainder of the 
image could then be subjected to the "blind" moisture 
classification algorithm. However, this approach cannot be 
expected to yield consistent results since for very dry soil 
moisture conditions many agricultural targets can appear 
similar to water (Figure 8a) or the water may be roughened by 
wind. In addition, for very wet soil conditions, many 
agricultural targets will be characterized by P r near nadir 
similar to that from the point targets (Figure 8c). A 
second, more satisfactory approach would be to incorporate 
si pr ior i knowledge of the spatial distribution of such 
features and filter them from moisture classification. This, 
of course, assumes the availability of a Level I land-use 
classification which could be scaled and rectified to the 
orbital radar imagery. 

In a similar fashion, the moisture classification error 
related to natural variability within the agricultural 
portions of the scene could be reduced if the radar data can 
be registered to topographic and crop distribution data. 

This would assume a mechanism for crop discrimination and 
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classification. In this case, each pixel element in the 
radar Image could be classified as to soil moisture using an 
algorithm tailored to be crop specific. Obviously, this 
approach is not currently feasible. 

However, since most of the confusion factors are 
spatially fixed and relatively invariant over short periods 
of time (excepting wind conditions), their effects on 
moisture classification accuracy can be minimized more 
economically by the multi-temporal change detection approach. 
In this technique, the radar imagery acquired at two dates 
are coregistered and their ratio yields a map of scene 
change. This process has been shown to be relatively simple 
to implement with L-band orbital imagery obtained by Seasat-A 
[8]. For a constant imaging geometry on the two dates (angle 
of incidence and azimuth view angle), the backscattered power 
received from cultural targets should remain approximately 
constant and that received from water bodies should remain 
nearly constant depending upon local wind conditions. Hence, 
these features should display little or no change in the 
multidate ratio images. On the other hand, all scene 
elements subject to change in backscatter category (such as 
planting, harvest, and tillage of agricultural fields) and/or 
subject to change in near-surface soil moisture status will 
yield a corresponding change in the multidate ratio images. 

If the time separation in multi-date observation is short 
relative to changes in crop development, then changes 
apparent in the ratio images will reflect relative moisture 


change and/or field status change related to tillage 
operations. Since surface slope is constant over the time 
interval, row direction is time constant in the absence of 
tillage, and surface roughness decays only slowly with time, 
the impact of these confusion factors upon the ratio of 
multidate received power should be negligible. 

The soil moisture distributions and the radar imagery 
simulated for Julian days 150 and 160 (wet and dry, 
respectively) are used to evaluate the utility of change 
detection for monitoring relative change in near-surfce soil 
moisture. The change in actual soil moisture between 

the two dates is shown in Figure 21. The graytone values in 
the image are linearly scaled to the difference function 
given as: 

AM fs = M fs (150) - M fs (160) (5) 

where the value in parentheses refers to Julian date. In 
producing Figure 21, a constant value of 128 (of 256 maximum) 
was added to AM fg , hence medium gray values such as those for 
the Kansas City area denote no change in soil moisture, 
bright areas denote considerable drying over the 10 day 
period, and dark areas denote an increase in near-surface 
(0-5 cm) soil moisture. In general, Figure 21 shows that 
drying conditions are prevalent over most of the test site 
except for scattered areas located primarily in the western 
portion (left side) due to rainfall (see Figure 5). 

Multidate registration of the radar imagery simulated at 
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increase in soil moisture. 




each of the three resolutions yields difference images which 

A 

are scaled to AM*, - via the blind classification algorithm 
t s 

(Equation 2). Image presentations of predicted change in 

A 

soil moisture AM r are shown for each radar resolution in 
f s 

Figure 22. In general, the direction (wetting or drying) and 
the magnitude of the true change in soil moisture observed in 
Figure 21 are faithfully reproduced for all radar 
resolutions. A noteable exception to this can be observed at 
the bottom center of each image in Figure 22. The black area 
denotes a predicted increase in soil moisture which is not 
observed in Figure 21. This discrepancy is the consequence 
of saturated to partially flooded soil conditions on Julian 
day 150 and moist conditions an Julian day 160 for this area. 
Hence, actual soil moisture has decreased while that 
predicted shows an increase since under flooded conditions 
the radar backscatter models generally yield low values 
of P r comparable to that from a water body. 

The area distributions of actual moisture change AM £g 

A 

and that predicted from the radar imagery AM fa are plotted in 

Figure 23. The sharp spike in the AM.. distribution at zero 

t s 

change is related to cultural features and water bodies. In 
general, it is apparent that the distribution of predicted 

A 

moisture change AM^ g as derived from the 100 m resolution 

radar most closely approximates the actual AM C distribution. 

f s 

The spatial averaging of the coarser radar resolutions causes 
them to be less sensitive to relatively large local change in 
AM F and thus the magnitude and extent of such changes tends 
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(b) 1 km by 1 km radar resolution. 
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Z Area 



J 


to be underpredicted. 

The actual and predicted change in soil moisture can be • 
compared on a grid cell basis by registration of the images 
in Figures 21 and 22. This procedure is, of course, subject 
to the registration errors discussed earlier for single date 
moisture classification due to changes in image geometry and 
position. For each pixel', the error in predicting relative 
moisture change can be defined as: 


AM 


AM - AM. 
t s £s 


( 6 ) 


The spatial distribution of E.„ is shown for each radar 
resolution in Figure 24. The brightest area on the scale bar 
denotes regions where the absolute magnitude of E^ M is within 
+/- 10% of AM^ and as graytone decreases the areas 
correspond to lE^I limits of +/- 20%, +/- 30%, and +/- 40% 
respectively as shown on the scale bar. For the 100 m 
resolution radar, 30% of the area is correctly classified to 
within +/- 20% of AM^ q and greater than 90% of the area to 
within +/- 30% of AM f g . In addition, most of the residual 
error is randomly distributed except for some classification 
error of large magnitude which is related to offsets in 
mechanical image registration as exemplified by linear 
features such as roads. For degraded radar resolutions of 1 
km and 3 km, the magnitude of classification errors increase 
and are spatially associated with edges between backscatter 
categories . 

The comparative error in moisture-change extimates E^ M 
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Figure 24. Spatial distribution of difference between actual change 
in soil moisture and that predicted from multidate radar 
observation. 
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(b) 1 km by 1 km radar resolution 


(c) 3 km by 3 km radar resolution 









for the three radar resolutions is shown in Figure 25 for all 

1.34 million grid cell comparisons within the test site. The 

corresponding percent of total area (124 km by 108 km) with 

absolute classification error less than a given magnitude is 

plotted in Figures 26a. Obviously, the 100 m resolution 

radar exhibits superior classification accuracy. However, if 

only the moisture variant pixels are compared (excludes 

cultural features, water bodies, and woodland) the 

distinction between resolutions shown in Figure 26b is not 

statistically significant; 78% and 89% of the area is 

correctly classified to within +/- 20% and +/- 30% of AM- , 

rs 

respectively. 

The effect of geographic subregion on the above results 
is shown in Figure 27. For the 100 m resolution radar, the 
change detection analysis results in superior classification 
accuracies for areas characterized by gentle topographic 
relief (cropland and urban subregions). For the coarser 
radar resolutions shown in Figure 27b and a, two effects are 
noted. First the influence of edges related to variance in 
the magnitude of AM fg between adjacent backscatter categories 
causes classification accuracy for all subregions to decrease 
relative to that for the rangeland/pasture subregion which is 
characterized by large field sizes. Secondly, the absolute 
classification accuracy decreases as a function of resolution 
for all subregions except rangeland/pasture. The large 
relative field size within the rangeland/pasture subregion 
and the large percent area occupied by range and pasture 


AM f s-AM fs 


Figure 25, Magnitude of error in estimates of relative soil moisture 
change as a function of test'site area. 
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(61.5%) ia largely responsible for the increase in 
class if teat ion accuracy using 1 Jem by 1 km radar data 
relative to that obtained using 100 m by 100 m radar data. 

Representative values of classification accuracy within 
each subregion for an error magnitude of +/- 20% of are 

shown in Table 6. These values show that 73% to 83% of the 
area within any subregion can be correctly classified as to 
within +/- 20% of actual soil moisture change for 100 m by 
100 m resolution radar imagery. In addition, these values 
are generally superior to those obtained for single date 
moisture classification shown in Table 5. 

5.0 CONCLUSIONS 

This simulation study reconfirms prior results that 
relatively high single-date moisture-classification 
accuracies can be achieved from orbital radar operating at 
4.75 GHz with HH polarization and at incidence angles of 7° 
to 17° relative to nadir. Furthermore, this study shows that 
classification accuracy is optimized for radar resolutions 
smaller than the expected field-size distribution of extended 
targets; a nominal sensor resolution on the order of 100 m by 
100 m is found to yield the most robust classification 
results for the majority of tested conditions. In addition, 
prior results have been extended to show that expected 
moisture-classification accuracy for a given sensor 
resolution is not independent of general soil moisture 


li 

'1 

i . 
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TABLE 6. Percent Area Correctly Classified to Within +/- 20% 
of the True Change in Soil Moisture AMf a from 
Julian Day 150 to Julian Day 160 


Subregions 

All 

LOO m 

Pixels 
1 km 

3 km 

Moisture Dep. 
100 km l km 

Pixels 
3 km 

Cropland 

82.4 

71.4 

70.8 

83.3 

74.0 

73.6 

Urban 

83.4 

73.9 

70.5 

80.4 

69.8 

64.2 

Rangeland/Pasture 

73.7 

81.77 

79.1 

74.1 

84.8 

82.3 

Woodlan ! 

74.7 

60.7 

52.6 

72.8 

73.5 

73.5 

Full Scene 

78.3 

74.6 

72.0 

78.3 

78.3 

76.8 
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condition or of the geographical mix of land-use, field-size 
distribution, and local topography. Finally, the use of 
multi-date radar imagery to estimate relative change in 
near-surface soil moisture status is shown to substantially 
reduce classification errors related to the presence of 
cultural features and water bodies, the presence of variable 
crop-canopy covers, and local variability in topographic 
relief . 

Based upon this study, a reasonable approach for the 
purposes of soil-moisture sensing would be to obtain the data 
at a sensor resolution on the order of 100 m (with a large 
number of independent looks) and then degrade the resolution 
where necessary by post-detection processing to average the 
moisture classification errors associated with local slope in 
regions of variable topographic relief. In addition, 
multi-temporal change-detection analyses could also minimize 
classification errors controlled by topographic relief as 
well as those errors that are related to intra- and 
inter-crop variance in radar backscatter ing [8]. 

6.0 ACKNOWLEDGMENTS 


This research was supported by the National Aeronautics 
and Space Administration, Johnson Space Flight Center, 
Houston, Texas, Grant NCC 9-6. 


83 


« 


REFERENCES 




[1] Ulaby, F. T. , M. C. Dobaon, J. Stiles, R. K. Moore, and 
J. Hoitzman, 1982, n A Simulation Stud/ of Soil Moisture 
Estimation by a Space SAR," Photo. Eng. Rem. Sena. , 
Vol. 48, No. 4, pp. 645-660. 

[2] Dobaon, M. C., F. T. Ulaby, and S. Moe 2 zi, 1982, 
"Assessment of Radar Resolution Requirements for Soil 
Moisture Estimation from Simulated Satellite Imagery," 
RSL Technical Report 551-2, University of Kansas Center 
for Research, Inc., Lawrence, Kansas. 

13] Anderson, J. R. , et al., 1976, "A Land Use and Land 

Cover Classification System for Use with Remote Sensor 
Data," Geol. Surv. Prof. Paper 964, U.S. Government 
Printing Office, Washington, D.C. 

[4] Hodges, T., 1981, "U.S. Crop Calendars in Support of 
the Earlier Warning Project," LEMSCO-14674, Lockheed 
Eng. and Management Services Co., Inc. 

[5] Ulaby, P. T., M. C. Dobson, and D. R. Brunfeldt, 1983, 
"Improvement of Moisture Estimation Accuracy of 
Vegetation-Covered Soil by Combined Active/Passive 
Microwave Remote Sensing," IEEE Trans. Geosc. and 
Rem. Sens . , accepted for publication 1983. 

[6] Ulaby, F. T., M. C. Dobson, J. Stiles, R. K. Moore, and 
J. Hoitzman, 1981, "Evaluation of the Soil Moisture 
Prediction Accuracy of a Space Radar Using Simulation 
Techniques," RSL Technical Report 429-1, University of 
Kansas Center for Research, Inc., Lawrence, Kansas. 

[7] Hoitzman, J. C. , V. s. Prost, J. L. Abbott, and 
V. H. Kaupp, 1978, "Radar Image Simulation," IEEE 
Trans. Geo. Elect. , Vol. GE-16, No. 5, pp. 296-303. 

[8] Brisco, B., F. T. Ulaby, and M. C. Dobson, 1983, 
"Spaceborne SAR Data for Land-Cover Classification and 
"'hange Detection," Digest of the 3rd IEEE 

^ntl. Geoscience and Remote Sensing Symp. , San 
Francisco, CA, August 31 - September 2. 

[9] Kansas State Board of Agriculture, 62nd Annual Report 
and Farm Facta , Topeka, Kansas, 1979. 

[10] Missouri Farm Fact3, Missouri Crop ana Livestock 
Reporting Service, Columbia, Missouri, April 1982. 

[11] Peck, A. J., R. J. Luxmoore, and J. L. Stolzy, 1977, 
"Effects of Spatial Variability of Soil Hydraulic 
Properties in Water Budget Modeling," Water Resources 


84 


Reaearch , Vol. 13, No. 2, pp. 348-354. 


[12] Peak, A. J., R. McQuivey, T. Keefer, E. R. Johnoon, and 
J. Erekaon, 1981, "Review of Hydrologic Models for 
Evaluating Use of Remote Sensing Capabilities,” NASA CR 
166674, Hydrex Corp., Pairfax, Virginia. 

( 13 J Hildreth, W. W., 1981, "Comparison of the 

Characteristics of Soil Water Profile Models," NASA - 
15700, Lockheed Eng. and Management Services Co., Inc. 

[14] Jackson, R. D., 1972, "On the Calculation of Hydraulic 
Conductivity," Soil Sci. Soc. Amer. Proc. Vol. 36, pp. 

380-382. 

[15] Armstrong, B. L. , 1978, "Derivation of Initial Soil 
Moisture Accounting Parameters from Soil Properties for 
the National Weather Service River Forecast System," 

NOA Technical Memorandum NWS HYDRO 37, National Weather 
Service. 

[16] Eagleman, J. R., 1971, "An Experimentally Derived Model 
for Actual Evapotranspirat ion, " Aqr ■ Met. , Vol. 8, 

pp. 385-394. 

[17] Clapp, R. B., and G. M. Hornberger , 1978, "Empirical 
Equations for Some Soil Hydraulic Properties," Water 
Res . , Vol. 14, No. 4, pp. 601-605. 

[18] Doorenbos, J., and W. 0. Pruitt, 1977, "Guidelines for 
Predicting Crop Water Requirements," FA0 Irrigation and 
Drainage Paper 24, United Nations Pood and Agriculture 
Organization. 

[19] United States Department of Agriculture, 1970, 
"Irrigation Water Requirements," Technical Release 21, 
USDA/Soil Conservation Service. 

[20] Hershfield, D. M. , 1961, "Rainfall Frequency Atlas of 
the U.S.," Hydrologic Service Division, Soil 
Conservation Serv'-e, U.S. Department of Agriculture, 
Washington, D.C. 

[21] Dunin, F. X., 1976, "Its Simulation for Field 

Conditions," Chapter 8 in J. C. Rodda (ed.): Facts of 

Hydrology , pp. 199-227, Wiley, London. 

[22] Komp, E. E., V. S. Frost, and J. C. Holtzman, 1983, 
"Reference Manual for the Radar Image Simulator," RSL 
Technical Report 581-2, University of Kansas Center for 
Research, i’nc., Lawrence, Kansas. 



APPENDIX A 


DYNAMIC SOIL WATER ACCOUNTING MODE L 

The purpose of a soil water-budget model within the 
context of realistic radar image simulation is to generate a 
distribution of near-surface (0-5 cm) soil moisture 
conditions at the spatial scale of the static terrain data 
base (100 m x 100 m) which responds to both static conditions 
(soil type, cover type, and surface slope) and dynamic 
conditions (crop stage, rain, and potential evaporation) on a 
time scale relevant to both the dynamics of the process and 
the orbital mechanics of an imaging satellite (daily basis). 
While many excellent water-budget modelg are available for 
various applications in agronomy and hydffology [11 to 15], no 
single model meets all the above criteria. Indeed, most such 
models require more detailed information on soil profile 
characteristics and weather conditions than is readily 
available for the simulation area. In addition, most models 
are designed to operate at a spatial scale much less than 
field size and over time increments significantly less than 
one day, or conversely, they are moat appropriately applied 
to very coarse integration times on the order of weeks for a 
3imple set of input parameters and at a macroscopic level 
much larger than field size. 

Because of the large size of the data base 
(approximately 1,339,000 grid cells), it is necessary to 
tailor a model that emphasizes the surface horizon and 
requires a minimum of information as to soil profile and 
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detailed local weather conditions/ and yet is still sensitive 
to daily variation in soil moisture. A schematic of the 
final process model is shown in Figure A.l; it consists 
largely of the following components: 
storm model, 
surface runoff model, 
crop development submodel, 
evapotranspiration model, and 
and interlayer redistribution model. 

When given dynamic inputs of crop type, crop stage of 
development, rainfall, and potential evaporation, the model 
acts upon the static terrain model to yield daily projections 
of 0-5 cm soil moisture for each grid cell. It also governs 
the redefinition of canopy cover categories based on crop 
calendar changes or local flooding conditions, and these 
categories are then used as input to the radar simulation 
program's target/sensor interaction model. 

A. 1 Storm Model 

Daily rainfall measurements a? reported by 25 stations 
located in and around tiirf test site were used as the basis 
for the storm model. Figure A. 2 shows the location of the 
test site. Table A.l shows the daily rainfall reported at 
each of these stations for the simulation period; May 18 
through June 9. A grid map of estimated rainfall, with a 
resolution of 3 km by 3 km, was produced from measured 
rainfall data at these irregularly spaced recording stations 



Figure A.!.. Dynamic Soil Water Accounting Model (SWAM). 
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for over/ day during the simulation period that all or part 
of the test site received some rain. Figures A. 3 and A. 4 
shows the amount of rain reported by each station on Julian 
day 144, and the estimated rainfall map for that day , 
respectively. These generated rainfall grid maps made 
available the total daily rainfall in cm received by each 
test site data base cell. An image representation of all 
rainfall grid maps has been shown on Figure 5. 

Rainfall intensity is calculated as a daily constant 
from the minimum recorded daily storm duration according to 


I o 10 a t * d. bt 
day day 


(A.l) 


where 


*day = dail y constant intensity, cm/hr 

daily minimum recorded duration, hrs. 
t = storm type (2-year or 5-year), and 

a and b are constant for each storm type. 

The constants a^. and b fc are solved from a plot of local 
rainfall intensity-vs-duration curves for recurrence 
intervals of 2 or 5 years. For each day of the simulation, a 
rainfall event is classified as either a 2-year or a 5-year 
event based upon the maximum recorded rainfall at all gauging 
stations on that day. If net daily rainfall at any gauge 
exceeds a critical value M, then that day will be classed as 
a 5-year event and a fc and b fc will be used from the 5-year 
inten3ity-vs-duration curve; otherwise a t and b fc will be used 
for a 2-year event. M is defined by 

M = 10 a * D b+1 /(b+l) (A. 2) 



where a and b a;;'e 2-year coefficients. Por the rainfall data 
given in Table A.l, the maximum net daily rainfall never 
exceeded M r therefore the 2-year coefficients were used in 
all precipitation events. 

A. 2 Surface Runoff Model 

The surface runoff model considers only the net effect 
of local surface slope and does not explicitly account for 
water retention and impoundment by soil surface roughness, 
tillage practices, and the presence of terraces. The water 
available for drainage as lateral surface flow is equal to 
the sum of standing water remaining from the previous daily 
accounting period plus the incident rainfall in excess of 
that which can infiltrate the surface layer and the root 
layer. The drainage D is computed from remaining standing 
water and local surface slope by 

D =» SW * (1.1 - 0 . 8 a ) (A. 3) 

where 

SW = standing water 

a ■=• the slope angle of the surface from horizontal in 
degrees. 

The term 1.1 - 0.8** is defined as the drainage coefficient 
and is plotted versus surface slope (in percent) in Figure 
A. 5. 


A. 3 Evapotranspiration Model 

Evapotranspiration is calculated differently for cropped 






and bare soil surfaces. For bare soil surfaces, the actual 
evaporation is depleted solely from the soil surface layer, 
while for vegetated surfaces a static root distribution model 
removes 30 percent of the actual evapotranspirat ion from the 
0-5-cm layer and removes the remaining 70 percent of actual 
evapotranspirat ion from the "root zone." For simplicity, the 
"root zone" is assumed to be one meter in depth and is 
treated as a constant with time and for all crops. 

For bare soil, actual evaporation, A E, is computed from 
potential evaporation, PE, as limited by antecedent soil 
moisture in the surface layer and soil hydraulic properties. 
Accounting is performed on a daily basis using the mean daily 
pan evaporation recorded at 11 stations in the study area as 
shown in Table A.l for 1981. 

An experimental model is used to calculate actual 
evaporation from potential evaporation PE: 

^ " PE * k soil * k storm < A ' 4 > 

where 


storm 


soil 


T 

PE 


lc_ 


- (24 - T)/24 , 

« soil limiting coefficient 
=> the duration of storm, and 

j * E pan' 


(A. 5) 


(A. 6) 


where 


k = pan coefficient, and 
P 

E = measured pan evaporation, 
pan 

The soil limiting coefficient k 0oil is defined by an 
experimental model [16] dependent upon PE and soil 
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properties . 


(A. 7) 


^soil * A + B <M R ) + C(MR) 2 + D(MR) 3 
where A, B, C, and D are empirically derived coefficients 
dependent upon PE, and MR is the moisture ratio. Regression 


fits to experimental data yield [16}: 

A - -0.05 + 0.732/PE (A. 8) 

B * 4.97 - 0.661 PE (A. 9) 

C - -8.57 + 1.56 PE (A. 10) 

D - 4.35 - 0.88 PE (A. 11) 

The moisture ratio MR is related to soil water retention 
characteristics via 

MR - (B - WP)/(FC - WP) (A. 12) 

where 


8 = measured soil moisture, 

WP => soil moisture at wilting point, and 
PC ~ soil moisture at field capacity. 

Assuming wilting point and field capacity to be defined as 
matric potentials of 15 bars and 1/3 bars, respectively, WP 
and PC can be defined from soil textural components using the 
approach of Clapp and Hornberger [17] 

FC " W 333 > Vb ' and (A. 13) 

WP = P (r Q /15,000) 1/b (A. 14 ) 

9 

where 


8 * soil moisture at saturation, 

9 

y » matric potential at saturation, and 

3 

b » an empirically derived value related to soil 
texture. 
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Por a given soil, 0 is calculated from the soil bulk density 

o 

profile and f and b are defined by A-horizon soil texture 
s 

using values given in [17]. Thus, for a given day, the terms 
in Eq. A. 4 are dependent on the antecedent soil moisture and 
the gross water-retention character i sties of each soil. 

For vegetated soil, the actual evapotranspiration, 
ET cro p, is computed by a modification of the Blaney-Cr iddle 
formulation used in estimating crop irrigation requirements 
[18,19]. Although the method is designed for an effective 
integration period of weeks to months, the simplicity of its 
input requirements makes this a practical approach for such a 
large number of coarse grid cells. Basically, crop 
consumption of water over the rooting depth varies with 
temperature, length of day, available soil moisture, crop 
type, crop stage of growth, relative humidity, and windspeed. 
To simplify the formulation, average measured values of 
temperature, day length, relative humidity, and windspeed are 
assumed on a seasonal basis for the simulation area. The 


resultant expression for ET becomes: 

crop 


? “PE * k * k , , 

crop crop storm 


(A. 15) 


where 


k^_^ = crop coefficient, 

crop c 

Crop coefficient as adjusted for mean local climate is 
plotted in Figure A. 6 as a function of number of days after 
planting for several of the crop covers included in the data 
base. Crop consumption of the water is seen to be dependent 
on both crop and stage of crop development. Before the crop 
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canopy has attained 20% ground coverage and again after 
harvest, the soil is treated as bare for both 
evapotranspirat ion and also for radar backscatter category. 

A. 4 Crop Development Model 

The length of time required for a given agricultural 
field in the simulation data base to progress from one 
crop-development stage to the next is established from data 
gathered by the Statistical Reporting Service of the United 
States Department of Agriculture. The simulation area lies 
at the East Central reporting district of Kansas (No. 6). 
Figure A. 7 presents a summary of mean crop development over a 
10-year period as enumerated by AgRISTARS [20] for this crop 
reporting district. These percentages are used to define 
crop development 3tage within the simulation on a 
field-by-field basis. Thus, each distinct agricultural field 
in the data base is assigned one of the 10 planting dates. 
Hence, there are ten different absolute crop calendars 
possible for each crop type identified in Table 1. 

Planting dates are randomly assigned to field codes for 
a specific crop based upon Figure A. 7. This procedure 
results in the introduction of a significant source of 
between-f ield variance of soil moisture within a given crop 
type due to the effect of crop development stage on 
evapotranspiration. It also allows for a given field to have 
its taget classification changed in Table 1, since a 
medium-rough bare field becomes a cropped field after 
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emergence, and finally reverts to bare soil status after 
harvest. As implemented, this procedure gives the data base 
a dynamic crop-category mix that can be modified to match 
regional agricultural practices such as double-cropping or 
dynamic soil surface roughness conditions. 

A - 5 Interlayer Water Redistribution 

Infiltration of water into the surface layer, 
percolation of water into the root zone, and capillary 
recharge of eurface layer moisture are controlled by the 
matr ic-potential profile as limited by soil structure. 

A pixel’s infiltration capacity during rainfall is given 
by [21J 

i t = t 1/2 f2 k s ir f ( 0 B - 6)) 1/Z + k s Z/2 (A. 16) 

where 

t = duration of rain event, hrs. 
k g - hydraulic conductivity at saturation 
t B ® suction at field capacity 
fl =■ proaity = 1 - p./p a 

® soil bulk density, g/cm 

3 

p g = soil specific density ® 2.65 g/cm for all soils. 
After rainfall ceases, infiltration proceeds at a rate 
defined by k_/2 for the remaining time of the accounting 
period (24-t) or until all standing water is depleted. Thus, 
a pixel's infiltration capacity from standing water is 
defined as 

i 8W - k a /2 * (24 - t) (A. 17) 
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and la limited by the amount of standing water. Hence, total 
infiltration into the surface layer of the soil, i fc , is 
determined by 

" *r + l ew < A * 18 > 

where 

i f < total rainfall received by the pixel 
i . , < standing water available. 

Water will percolate from the surface layer (0-5-cm) 
into the root zone for all accounting periods where the 
surface layer’s water content after infiltation exceeds the 
water content at field capacity (as determined by Equation 
A. 13), such that final surface-layer's water content is 
reduced to less than porosity. 

This is accomplished by first allowing excess water to 
drain from the root zone (5-100 cm depth). One third of the 
volumetric moisture in excess of root zone field capacity is 
allowed to drain gravitationally each day and hence is 
removed from further accounting periods. Then, assuming that 
the water content in the surface layer exceeds field 
capacity, the excess is permitted to .percolate into the root 
zone at the minimum of either 

T k 

Rj •= (A. 19) 

1 2a 

where 

*» net percolation into the root zone, 
a ■= a damping coefficient arbitrarily set to 48, and 


T - duration of accounting period - 24 hours 
or 1/3 of excesB water is allowed to percolate 

R 4 - (6 - FC)/3 (A. 20) 

where B and FC are for the surface layer. 

When evapotranspirative losses cause surface-layer water 
content to be reduced below wilting point, capillary recharge 
of the surface 5 cm of soil is allowed to occur during the 
night for a duration of 12 hours. The rate of the surface 
recharge is equal to k g /2 and is arbitrarily limited to a 
maximum of 0.25 cm of water. Furthermore, capillary recharge 
is not allowed to raise surface layer water content above 
wilting point. 

A . 6 Within-Field Variability in Surface Soil Moisture 

Prior to radar image simulation, the surface layer soil 
moisture values determined by the water-budget model for each 
100-m by 100-m grid cell are randomized to approximate the 
natural variability in soil moisture measured within 
"homogeneous" fields. Randomization was performed on a 
grid-cell basis by a Gaussian random-number generator with a 
standard deviation of 6 percent M fg [2]. 

A. 7 Generation of Soil Moisture Distributions 

The dynamic soil water accounting model (SWAM) was 
initialized on Julian day 138 and moisture distribution maps 
of the test site were produced for every day of the 
simulation period. 
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These moisture maps indicated the percent of the 1/3-bar 
water content M^ s in the 0-5 cm layer where 

M fg - 100 x 0/PC (A. 21) 

where 

e - measured soil moisture 

PC = soil moisture at field capacity. 

The resultant distributions were then examined and the three 
most closely approximating moderately dry, moist, and wet 
soil surface conditions were selected for radar image 
simulation. Image representation of 0-5-cm soil moisture 
distribution for Julian day 141, 150 and 160 are shown on 
Figure 6. 


DEFINITIONS 


PAUSDATE 

STRTDATE 
STOPDATE 
RA INDATE 

ALL I NTS 

RAIN 

INTNSITY 

DUR 

PERCENTS 

SFBULK 

RTBULX 

B 

FSUCTION 

SSUCTION 

SHYDCOND 


An array which contains the Julian da/3 
on which the output moisture map needs 
to be saved. 

Julian day on which the process should 
begin. 

Julian day on which the process should 
stop. 

An array containing the Julian dates which 
all or part of the database received some 
rain. 

An array containing the mean rain intensity 
of each rainy day. 

Amount of rain received by a cell on a 
certain day in cm. 

Intensity of the rain for a cell in cm/hour. 

Duration of the rain for a cell in hoot 3 

In soil data subroutine. An array of 
percent probability of occurrence of soil 
bulk density associated with each of 
eight soil types present in our data base. 

Quantized levels of surface layer (0-5 cm) 
bulk density associated with "PERCENTS". 

Quantized levels of root layer (5-100 cm) bulk 
density associated with "PERCENTS". 

An array containing b values for all 
15 soil textures as estimated from Clapp & 
Hornberger, 1978. 

An array containing the suction (at field 
capacity) for all 15 soil textures (see Clapp & 
Hornberger, 1978). 

An array containing the suction (at 
saturation) for all 15 soil textures. 

An array containing the hydraulic 
conductivity at saturation k $ for all 15 soil 
textures. 
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SATHC 

PDATES 

STAG EDA ¥ 

KEQCONST 

KCROP 

SW 

SWINF 

RAININF 

MFC 

DRAIN 

SWRUNOFF 

RECHRG 

ETO 

SWEVAP 

KSOIL 


Hydraulic conductivity at saturation. 

An array containing ton different planting 
dates for each crop type. 

An array containing the number of days 
after the planting date which the crop 
advances to a new crop growth stage (five 
different stages) for each crop type. 

An array containing two parameters 
(slope and intercept) describing the 
change in K crop at each stage and for 
each crop type. 

Crop transpiration coefficient. 

Standing water (cm). 

Amount of standing water which 
infiltrates to the surface layer (cm). 

Amount of rain which infiltrates 
to the surface layer. 

Water content expressed as a percent 
of field capacity. 

Amount of excess water which is drained 
from the root zone (cm/cm). 

Amount of water runoff from standing 
water (cm) . 

Capillary recharge (cm). 

An array containing the potential 
evaporation (cm) for every day of the 
simulation period. 

Amount of evaporation from standing ' 
water . 

Bare soils evaporation coefficient. 
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T h l f °H 0Win 9 Var iablea are prefixed by "SF" or "rt" ' 

«mt S»pcctW„l^ f£,Ca l£y " <0 ' 5 cm > «ne (5-100 

soil's bulk density 
soil's porosity 

soil's water content at field capacity (cm/cm) 
soil's water content at wilting point (cm/cm) 
water content (cm/cm) 
amount of evaporation (cm) 


BO 

PROS 

FC 

WP 

WC 

EVAP 
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UNIVERSITY OF KANSAS REMOTE SENSING LAB 

PROCRAM SUITE 

l RADAR 

SIMULATION 

REF. « t RSL 

REPORT 601-1 

PROGRAM NAPIEt 

language i 

SWAM 

FORTRAN 

AUTHOR l 
77 

SAIED MOEZZI 

DATE* MAY 1933 


PURPOSE t THE PURPOSE OP PROGRAM SUAM (SOIL UftTER ACCOUNTING MODEL) 
IS TO GENERATE A DI5RIBUTI0N OF NEAR-SURFACE (0-5 CN> SOIL MOISTURE 
CONDITION AT THE SPATIAL SCALE OF THE STATIC TERRAIN DATA BASE WHICH 
RESPONDS TO BOTH STATIC CONDITIONS (SOIL TYPE, COVER TYPE, AND SURFACE 
SLOPE) AND DYNAMIC CONDITIONS (CROP STAGE, RAIN, AND POTENTIAL EVAPO- 
RATION) ON A DAILY BASIS. 


NAME 

IROUl 

IR0U2 

ICOU 

ICOL2 

IOTCOL 

NCOL 

NPAUSE 


NAME 

RAINFALL 


INTRLAYR 


INITIALZ 

EVAPORAT 

CELLDATA 

CROPDATA 

SOILDATA 

GETFILES 

IOCALLS 


UPDTHIST 

OTPTHIST 


PARAMETER DEFINITION 

DESCRIPTION 

FIRST ROU OF THE INPUT MATRIX TO BE PROCESSED 
LAST ROU OF THE INPUT MATRIX TO BE PROCESSED 
FIRST COLUMN CELL OF THE INPUT ROUS TO BE PROCESSED 
LAST COLUMN CELL OF THE INPUT ROUS TO BE PROCESSED 
NUMBER OF CELLS IN EACH OUTPUT ROW 
NUMBER OF CELLS IN EACH INPUT ROU 
NUMBER OF TIMES THAT PROGRAM SHOULD PAUSE DURING 
SIMULATION PERIOD FOR SAVING THE MOISTURE MAP 


I 


SUBROUTINES REQUIRED 


DESCRIPTION 


56 COPNFIL 

57 C 


IRETURNS THE AMOUNT OF RAIN (CM), DURATION (HOURS) AND 
(INTENSITY < CM/HOUR > FOR A GIVEN CELL ON A SPECIFIED 
(JULIAN DAY. 

(THIS ROUTINE CONTAINS ALL 4 UATER ACCOUNTING MODELS 
(EACH AS A SEARATE ENTRY. THESE ENTRIES ARE l SURFINF, 
1 ROOTINF, RUNOFF AND RECHARGE. 

(THIS ROUTINE IS USED FOR INITIALIZATION PROCESS AND 
(HAS TUO ENTRIES. THESE ARE * COMMENCE AND DAUN. 

(THIS ROUTINE SIMULATES THE EVAPORATION PROCESS FOR A 
(GIVEN CELL ON A SPECIFIED JULIAN DAY. 

(THIS ROUTINE GETS ALL STATIC CONDITIONS OF A GIVEN 
I DATA BASE CELL. 

(THIS ROUTINE GETS THE DYNAMIC CONDITIONS OF A GIVEN 
t DATA BASE CELL UH1CH IS REGISTERED AS A CROP TYPE. 
(THIS ROUTINE GET ALL THE REQUIRED INFORMATION 
(THAT IS BASED ON THE SOIL TYPE FOR A GIVEN CELL. 

(THIS ROUTINE OPENS ALL THE INPUT AN OUTPUT FILES, 
(THIS ROUTINE CONTAINS TUO ENTRIES FOR READING AND 
I WRITING INPUT AND OUTPUT RECORDS. THESE ARE* READREC, 
(AND URITDATA. 

(THIS ROUTINE IS USED FOR UPDATING ft GIVEN HISTIGRAM. 
(THIS ROUTINE IS USED FOR WRITING OUT A GIVEN 
(HISTOGRAM. 

(EXTERNAL FUNCTION CALLED BEY 'GETFILE5 ' ROUTINE 
I 


58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 


PARAMETER ( IROUl ■ 1 , IROU2-1077, ICOL1-1 , IC0L2- 1245, IOTCOL- 1245 ) 
PARAMETER ( NPAUSE *4 , NCOL- 1245 ) 

INTEGER PAUSDATE( NPAUSE ), FC( 14 ), HHMMSS< 3 ) 

INTEGER WATER, SOIL, ELEU, CATG, CAT, COL, ROU, DATE 

INTEGER STRTDATE, STOPDATE. CROP* FIELD, DAY 

REAL KCROP, INTNSITY, MFC 

CHARACTER*8 TYPE, TYPENOU 


LOGICAL 


PAUSE 
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v'-zi ra 

OF P/Ou ("iAUTY 


?1 COMMON /BUF3/ MFCOUTl NPAUSE. NCOL >, ICATOUYt NPAUSE. NCOL ) 

72 SPECIAL COMMON BUF3 

73 COMMON /BUF4/ MFCHISTt NPAUSE , 256 J. ICATHIST< NPAUSE . 31 ) 

74 SPECIAL COMMON BUF4 

75 COMMON /BLOCK/ 

76 & SATHC, SUCTION, SFPROS, RTPROS, SFFC. RTFC, SFUP, 

77 l SU, SFUC, RTUC, KCROP, 

78 & DAY, RAIN, BUR, INTNSITV, SLOPE, TVPENOU 

79 COMMON /FILCOD/ FC 

80 DATA UATER/ 10 /. STRTDATE/ 138 /. STOPDATE/ ISO / 

8« DATA PAUSDATE/ 140, 141, 1S0, 160 / 

82 C 

83 C 
34 C 

85 URITEt 13,101 ) IROUt, IRQU2, ICOL1, ICOL2 

86 C 

87 C INITIALIZE ALL NECESSARV UARIABIES 

88 C 

89 CALL COMMENCE 

90 C 

91 C PROCESS EUERV CELL IN THE DATA BASE 

92 C 

93 FOR ROD* IROUl , IRCU2 

94 C 

95 C AFTER PROCESSING EVERY 100 RECORDS SEND A MESSACE TO TO TERMINAL 

96 C 

97 /,F< MODtROU, 10 ) ,EQ. 0 ) THEN 

98 CALL TIME(HHMMSS) 

99 URITE ( 11,103) ROU, HHMMSS 

100 END IF 

101 FOR COL* ICOL1 , 1C0L2 

tea c 

103 C GET REQUIRED INFORMATIONS FOR THE CELL BEING PROCESSED 

104 C 

105 CALL CELLDATA< ROU. COL. ELEU, SOIL, CATG. SLOPE, TYPE, CROP, FIELD > 

186 IF< TYPE .EG. 'NONAGRIC' ) THEN 

107 C 

108 C THIS IS NOT AN AGRICULTURAL CELL THEREFORE SHOULD NOT BE 

109 C TREATED IN MOISTURE COMPUTATION. UPDATE THE OUTPUT ROU AND 

110 C START UITH NEXT CELL IN THE DATA BASE. 

111 C 

112 FOR IP-1, NPAUSE 

113 NFCOUT ( IP , COL ) • 0 

114 ICATOUT< IP , COL > * CATG 

115 CALL UPDTHISTt 0 , MFCHIST, IP, 0, 250 ) 

116 ICAT * CATG / 10 

117 CALL UPDTHISTt ICAT, ICATHIST, IP, 0, 25 ) 

118 END FOR 

119 ELSE 
128 C 

121 C ELSE THIS CELL IS AN AGRICULTURAL TYPE, START THE MOISTURE 

122 C COMPUTATION AND CONTINUE FOR THE ENTIRE SIMULATION PERIOD. 

123 C 

124 C GET MORE INFORMATION ABOUT THE UNDERLAYING SOIL 

125 C 

126 CALL SOXLDATAt SOIL, SFPROS, RTPROS, SFFC, RTFC, SFUP, 

127 & SATHC, SUCTION ) 

128 C 

129 C INITIALIZE THIS CELL'S MOISTURE FOR DAY ZERO 

130 C 

131 CALL DAUN 

132 C 

133 FOR DATE- STRTDATE, STOPDATE 

134 C 

135 DAY • DAY ♦ 1 

136 C 

137 C IF IT IS A RAINY DAV, THEN GET AMOUNT, INTENSITY AND THE 

138 C DURATION OF THE RAINFALL ON THIS GROUND CELL BEING PROCESSED 

139 C 

140 CALL RAINFALLC DATE, ROU. COL, RAIN, INTNSITV, DUR ) 
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OF 


y.tVl 


; v- n 


141 C 

142 C 

143 C 

144 C 

145 

146 

147 

148 C 

149 C 

150 C 

151 C 

152 C 

153 C 

154 

155 C 

156 C 

157 C 

158 

159 C 

160 C 

161 C 

162 

163 C 

164 C 

165 C 

166 

167 C 

168 C 

169 C 

170 

171 C 

172 C 

173 C 

174 C 

176 

177 C 

178 C 

179 C 

180 C 

181 C 
122 C 

183 

184 

185 

186 

187 

188 

189 

190 100 

191 C 

192 C 

193 C 

194 C 

195 

196 C 

197 C 

198 C 

199 

200 C 

201 C 

202 C 

203 

204 C 

205 C 

206 C 

207 C 

208 C 

209 C 

210 C 

> 


IF THE GROUND CEU IS REGISTERED AS A CROP THEN GET 
KCROP AND CROP STAGE 

IF( TVPE .EG. 'CROP ' ) THEN 

CALL CROPDATA< CROP, FIELD, DATE, KCROP. TYPENOU > 

END IF 

PROCESS ALL SOIL UATER ACCOUNTING PIODELS 

It PONDING AND INFILTRATION INTO THE SURFACE LAVER (0 - 5 CPI) 
CALL SURFINFL 

2t PERCOLATION OF UATER INTO THE ROOT ZONE (5 - 95 CM) 

CALL ROOTINFL 

3 ■ STANDING UATER RUNOFF DUE TO LOCAL SLOPE 

i 

CALL RUNOFF 
4 1 EUAPOTRANSPIRATION 
CALL EUAPORAT 

5» CAPILLARY RECHARGE OF THE SURFACE LAYER 
CALL RECHARGE 


COMPUTE’ X OF FIELD CAPACITY OF SOIL flOISTURE BASED ON SURFACE 
LAVER'S UATER CONTENT 
NFC « 100.0 S SFUC / SFFC 


CHECK TO SEE IF THIS IS A PAUSE DAY. IF IT IS THEN RECORD 
THE COPIPUTED *NFC* AND THE REASSIGNED CATEGORY. 

(PAUSE DAY IS THE DAY THAT THE PIOISTURE NAP NUST BE SAUED) 

PAUSE - .FALSE. 

FOR IP-1, NPAUSE 

IF < DATE .EG. PAUSDATE( IP ) ) THEN 
PAUSE * .TRUE. 

GOTO 100 
END IF 
END FOR 

IF ( PAUSE ) THEN 

APPLY A GAUSIAN DISTRIBUTION UITH COMPUTED "MFC 1 AS THE 
PIEAN. AND 6X P1FC AS THE STANDARD DIUIATION 

IP1FC - NINT< RANNC P1FC , 6.0 ) ) 

SET LOUER LIMIT OF X FIELD CAPACITY TO ONE 

IF( IP1FC .LT. 1 ) IPIFC-1 

SAUE COMPUTED ’MFC* FOR THIS PAUSE DAY 

MFCOUTC IP , COL ) • IMFC 

REASSIGNMENT OF THE REGISTERED CELL'S CATEGORY 
It CHANGE THE CATEGORY TO UATER IF THERE IS STANDINNG UATER 
ON THIS GROUND CELL 

2t CHANGE THE CATEGORY TO UATER IF SURFACE LAVER'S UATER 
CONTENT EXCEEDS THE UNDERLAYING SOIL'S PROSITY 
3* IF THERE IS NO STANDING UATER AND CELL IS REGISTERD AS A 


no 


C >’ i . ... 

of pour Qj/vurr 


ail c 

CROP TYPE, THEM CHANGE THE CATEGORY TO DARE SOIL WHEN 

SIS c 

BEFORE EMERGENCE OF THE CROP OR AFTER HARVEST 

£13 C 



214 

IF( SU .GT . 0 ) THEN 


215 

CAT - HATER 


216 

ELSE 


217 

IF< SFUC .LT. SFPROS > THEN 


218 

CAT • CATG 


219 

IF < (TYPE. EG. 'CROP ') . AND. (TVPENOU.EQ. 'SMTHBARE' ) ) CAT • 210 

220 

IF( (TYPE. EO. "CROP ') .AND. (TYPENQU.EQ. 'MEDMBARE' ) ) CAT • 200 

221 

ELSE 


222 

CAT • HATER 


233 

END IF 


224 

END IF 


225 

ICATOUT< IP , COL ) • CAT 


226 

ICAT • CAT / 10 


227 C 



228 C 

UPDATE THE HISTOGRAMS 


229 C 



230 

CALL UPDTHISTl IMFC, MFCHIST, IP, 0, 250 ) 

t 

231 

CALL UPDTHIST( ICAT. ICATHIST, IP, 0. 85 > 


232 C 



233 

END IF 


234 

END FOR 


235 C 



236 C 

DONE UITH MOISTURE ESTIMATION FOR THIS CELL 


237 C 



238 

END IF 


239 

END FOR 


240 C 



241 C 

DONE UITH ALL THE COLUMNS OF THIS ROU 


242 C 

URITE OUT THE COMPUTED MFC AND THE REASSINED 

CATEGORIES OF THIS 

243 C 

ROU TO THE OUTPUT FILES FOR ALL PAUSE DATES. 


244 C 



245 

call uritdatac i. hfcout. iotcol j 


246 

CALL URITDATAt 2, ICATOUT, IOTCOL ) 


247 C 



248 

END FOR 


249 C 



250 C 

SOIL MOISTURE ESTIMATION IS DONE FOR THE ENTIRE DATA BASE 

251 C 

URITE OUT A REPORT OF THE FINAL MOISTURE AND 

CATEGORY HAPS 

252 C 



253 

FOR IP-1, NPAUSE 


254 

URITE ( 13, 104 ) 'MFC ", PAUSDATE( IP ) 


255 

CALL 0TPTHIST( MFCHIST , IP, 0, 250 ) 


256 

URITEt 13,104) 'CATEGORY', PAUSDATE* IP > 


257 

CALL 0TPTHIST1 ICATHIST, IP, 0, 25 ) 


258 

END FOR 


259 

CALL TIME(HHMMSS) 


260 

URITE ( 13, ' ( 1 COMPLETED AT *,3A3>') HHMMSS 


261 

URITEC 11, ' < * COMPLETED AT *,3A3>') HHMMSS 


262 

URITE( 11,102) 


263 

URITEt 13, 102 ) 


264 

STOP 


265 101 

FORMAT (/V/, ' SOIL UATER ACCOU 

N T I N G', 

266 

& ' P R 0 G R A M'//,' SHAM UAS PROCESSED ON 

THE DATA BASE'// 

267 

& ' ROU' 15, ' THROUGH ROU'IS, ' , COL'IS' 

THROUGH COL'IS) 

268 102 

FORMATt' * * * ALL DONE***') 


269 103 

FORMAT!' PROCESSED THROUGH 

R E C 0 R D ' , 5X, 14 

270 

& 3X,3A3 ) 


271 104 

FORMAT ( * 1 • , A8, ' HISTOGRAM FOR JULIAN DAV'I5> 

272 

END 
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273 

274 

275 

276 

277 

278 

279 

280 
281 
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284 

285 
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287 

288 

289 
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319 

320 

321 
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323 

324 
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326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

c 


10 


c 

99 


OF FOCrt 


THIS ROUTINE READS THE RAINFALL DATA AND RETURNS THE AMOUNT 
OF RAIN, DURATION, AND INTENSITY FOR A GIVEN CELL UITHIN THE 
DATA BASE ON A SPECIFIED JULIAN DAY. 

RAIN IS IN UNITS OF CENTIMETERS , INTENSITY IS IN CM/HOUR, 
AND DURATION IS IN HOURS. 


SUBROUTINE RAINFALL * DATE, DBROU, DBCOL, RAIN. INTNSITY, DUR ) 
IMPLICIT INTEGER ( A - 2 > 

PARAMETER ( NDAY-13 , NRCOL-42 ) 

COMMON /BUF5/ ROURAIN* NRCOL.NDAY), RAINDATE ( NDAV ) 

SPECIAL COMMON BUFS 
DIMENSION FC< 14 ) 

COMMON /FILCOD/ FC 

REAL RAIN, DUR, INTNSITY. ALLINTS* NDAY > 

DATA RECPTR / 0 / 

DATA RAINDATE 

& ✓ 138, 139, 143, 144, 148, 149, 150,151, 153, 154, 155, IS9, 161 / 

DATA ALLINTS 

& / 1.2,1. 7, 1.7, 3. 6. 2. 3, 2. 3, 3. 6. 1.4, 3. 6, 3. 6, 3. 6, 2. 3, 2. 3 / 


READ RAINFALL DATA FOR ALL RAINY DAYS FOR THIS GROUND CELL 

RAINROU • ( ( DBROU - 1 ) / 30 + 1 ) 

RAINCOL • < * DBCOL - 1 > / 30 4 1 ) 

UHILE ( RECPTR .LT. RAINROU > 

FOR COL-1, NRCOL 

READ* FC(5>, IOSTAT-IOS) * ROURAIMCOL, DAY ), DAY-1, NDAV J 
IF ( IOS .NE. 0 ) GOTO 99 
END FOR 

RECPTR • RECPTR + 1 
END UHILE 


CHECK IF THE DATE GIVEN UAS A RAINY DAV 
FOR DAY-1, NDAV 

IF* RAINDATE* DAY ) .EQ. DATE ) GOTO 10 
END FOR 

NO STORM ON THIS DAY, RETURN TO THE CALLING PROGRAM 

RAIN * 0. 

INTNSITY • 0, 

DUR • 0 . 

RETURN 

A STORM OCCURED ON THIS DAY, GET THE AMOUNT OF RAIN RECIEVED 
BY THIS CELL ON THE GIUEN JULIAN DAV 

RAIN • REAL* ROURAIN* RAINCOL, DAV ) >/ 10.0 
INTNSITV • ALLINTS* DAV ) 

DUR ■ RAIN / INTNSITV 
RETURN 

URITEC11,'*1X, 'ERROR ttxt UHILE READING RAINFALL ')') 

STOP 

END 
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400 C 

401 C 
408 C 

403 

404 
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407 C 
488 C 
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u! ; g'j.aJiV 


INTRLAYR 


THIS ROUTINE CONTAINS AU 4 UATER ACCOUNTING MODELS, EACH 
AS A SEPARATE ENTRY. THESE ARE 1-SURFACE INFILTRATION. 
a-ROOT INFILTARION, 3-RUNOFF. 4-RECHARCE 

SUBROUTINE INTRLAYR 

INTEGER FC< 14 ), DAY 

REAL KCROP, INTNSITY, flAXRTINF, MINRTINF 

CHARACTERI8 TVPENOU 
COMMON /BLOCK/ 

& SATHC. SUCTION. SFPROS. RTPROS, SFFC. RTFC , SFUP. 

& SU. SFUC . RTUC, KCROP, 

l DAY, RAIN, DUR, INTNSITY, SLOPE, TYPENOU 
COnnON /FILCOD/ FC 

ENTRY SURFINFL 


C0I1PUTE AMOUNT OF ‘RAIN* UHICH INFILTATES TO SURFACE LAYER 
IF( RAIN .GT. 0 > THEN 

Ri • SORT ( DUR > * SORT ( S*SATHC*SUCTION*( SFPROS-SFUC) ) ♦ 
t SATHC * DUR / 8 

RAININF • AMIN1 ( RAIN , RI ) 

SU ■ SU ♦ ( RAIN - RAININF ) 

ELSE 

RAININF • 0 
END IF 

CONFUTE AMOUNT OF 'STANDING UATER’ UHICH INFILTRATES TO SURFACE 
LAYER AFTER RAINFALL CEASES 

IF< SU .CT. 0 ) THEN 

SUI • SATHC « ( 84 - DUR > / 8 

SUINF • AHIN1 ( SU , SUI ) 

SU • SU - SUINF 
ELSE 

SUINF - 0 
END IF 

TOTAL AMOUNT OF INFILTRATION TO THE SURFACE LAYER IS THE SUN 
OF RAIN AND STANDING UATER INFILTARTION 

TOTALINF • RAININF + SUINF 

COMPUTE UATER CONTENT PER CENTIMETER OF SURFACE LAYER 

SFUC • SFUC + TOTALINF / 5.0 
RETURN 


ENTRY ROOTINFL 


COMPUTE .PERCOLATION OF UATER FROM THE SURFACE LAYER INTO 
THE ROOT ZOON ( 5-100 CM ) 

FIRST DRAIN THE EXCESS UATER OUT OF ROOT LAYER 

IF< RTUC .GT. RTFC ) THEN 
DRAIN - (RTUC - RTFC) / 3 
RTUC • RTUC - DRAIN 
END IF 

THEN PERCOLATE 
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rrv. 


409 C 

410 

41 1 
418 

413 

414 

415 

416 

417 

418 

419 

420 

421 

422 C 

423 C 

424 C 

425 C 

426 C 

427 
42B 

429 

430 

431 

432 C 

433 


'(3- C 

435 C 

436 

437 C 

438 C 

439 C 

440 C 

441 C 

442 


445 

446 

447 

448 

449 C 

450 C 

451 

452 C 

453 C 

454 C 

455 C 

456 C 

457 

458 

459 

460 

461 

462 

463 

464 

465 

466 

467 

468 

469 C 

470 

471 


IF< 5FUC .GT. SFFC ) THEN 
SFUC! • SFUC 
RTINFl • 0.25 * SATHC 
RTINF2 • ( SFUC - SFTC ) / 3 
SFUC • ANIN1 < SFUC! , SFPROS ) 

SFUC • ANAXK SFUC-RTINF2 , SFUC1-RTINF 1 ) 

RTINF • ( SFUC1 - SFUC) * 5 
RTUC » RTUC ♦ RTINF/95 
ELSE 

RTINF • 0 
RETURN 
END IF 

tlOU CHECK THE SURFACE UATER CONTENT, IF IT EXCEEDS THE 
PROSITV OF THE SOIt TVPE PtOUE THE EXCESS UATER TO THE 
STANDING UATER 

IF< SFUC .GT. SFPROS ) THEN 
SUADD • ( SFUC - SFPROS ) * 5 
SfUC • SFPROS 
SU * SU ♦ SUADD 
END IF 

RETURN 


ENTRV RUNOFF 


CONFUTE RUNOFF CAUSED BV THE SLOPE FOR ‘STANDING UATER* 


IF( SU .GT. 0 J THEN 
SURUNOFF » SU 8 l 1.1 - 


SU • ANIN1 ( SU 

SUR'JNOFF • 0. 
END IF 
RETURN 


SU-SURUNOFF > 


a.fl 115L0PE ) 


ENTRV RECHARGE 


CAPILLARV RECHARGE IS ALLOUED TO CcCURE DURING NIGHT FOR A 
DURATION OF 12 HOURS 

IF < SFUC .LT. SFUP ) THEN 

SFUCt • SFUC 

RECHRG • 0.25 

SFUC • SFUC1 + RECHRG/5 

I F < SFUC .GT. SFUP ) THEN 

SFUC ■ SFUP 

RECHRG • (SFUC - SFUC1) * 5 
END IF 

RTUC • ( 95 * RTUC - RECHRG ) r 95 
ELSE 

RECHRG - 0 
END IF 

' RETURN 
END 
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472 c 

473 C 

474 C 

475 C 

476 C 

477 C 

478 C 

479 C 

480 

481 C 

482 C 

483 

484 

485 

486 

487 

488 

489 

490 

491 C 

492 C 

493 

494 C 

495 C 

496 C 

497 C 

498 C 

499 C 

500 

501 

502 

503 

504 

505 

506 

507 

508 C 
569 C 

510 C 

511 

512 C 

513 C 

514 C 

515 C 

516 

517 

518 

519 

520 

521 C 

522 101 

523 
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THIS ROUTINE IS USED FOR INITIALIZING THE VARIABLES AS UELL 
AS INITIAL MESSAGES TO THE TERMINAL AND OUTPUT REPORT FILE. 

SUBROUTINE INITIALZ 


INTEGER DAYTIME! 2 >, DDMMMVVC 3 ), HHMMSSC 3 >, FCC14) 

REAL KCROP, INTNSITY 

CHARACTERS TYPENOU 
COMMON /'BLOCK/' 

& SATHC, SUCTION. SFPROS, RTPROS, SFFC, RTFC, SFUP, 

& SU, SFUC, RTUC. KCROP, 

& DAY, RAIN, DUR, INTNSITY, SLOPE, TYPENOU 
COMMON /FILCOD/ FC 


ENTRY COMMENCE 


AT THE BEGINING URITE OUT A MESSAGE TO THE TERMINAL AND 
GET A SEED FOR RANDOM NUMBER GENERATOR FUNCTIONS BASED 
ON THE COMPUTER CLOCK 

CALL DATE! DDMMMYY > 

CALL TIME( HHMMSS ) 

CALL JDATE! DAYTIME ) 

URITEUl.iei ) HHMMSS. DDMMMYY 
URITE! 13,101) HHMMSS, DDMMMYY 
PRIMNO • DAYTIME! 2 ) 

CALL IRANP( PRIMNO ) 

RETURN 


ENTRY DAWN 


THIS ENTRY INITIALIZES THE MOISTURE CONTENT OF A CELL 

SU - 0 
DAY • 0 
SFUC - SFFC 
RTUC • RTFC 
RETURN 

FORMAT!//' TIME! '3A3, ' , DATEt '3A3) 

END 


> 


DAY 
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sas c 

536 C EVAPORAT 

537 C 

538 C THIS ROUTINE SIMULATES THE EVAPOTRANSP ITAION ON A GIVEN 

539 C DAY BASED ON THE DATA GATHERED. BARE SOIL, CANOPY COVERED 

530 C OR UATER COVERED CELLS ARE EACH TREATED DIFFERENTLY. 

531 C 

533 SUBROUTINE EUAPCRAT 

533 C 

534 PARAMETER < NSDAY- 34 ) 

535 C 

536 INTEGER FCt 14 ), DAY 

537 REAL £T0( NSDAY ), MR, KSOIl, KSTORM, KSUEVAP, KCROP, INTNSITV 

538 REAL At NSDAY ), Bt NSDAY ), C( NSDAY ), D< NSDAY ) 

539 LOGICAL FRSTCALL 

540 CHARACTER*8 TYPENOU 

541 COMMON /DATA/ SU INF, RAIN INF, RT INF, SURUNOFF,RECHRG, SUE VAP, SFEVAP, 

543 & RTEVAP 

543 COMMON /BLOCK/ 

544 & SATHC, SUCTION, SFPROS, RTPROS, SFFC. RTFC, SFUP, 

545 & SU, SFUC, RTUC, KCROP, 

546 & DAY, RAIN, DUR, INTNSITY, SLOPE, TYPENOU 

54? connoN /filcoo/ fc 

548 DATA FRSTCALL/ .TRUE. / 

549 DATA ET0 / 0.48. 0.33, 0.31, 0.48, 0.61, 0.58, 0.40, 0,46, 

550 6 0.38, 0.37,-0.33, 0.15, 0.3S, 0.35, 0.43, 0.33, 

551 & 0.38, 0.38, 0.31, 0.31, 0,40, 0.68, 0.55, 0.59 / 

552 C 

553 IF ( FRSTCALL ) THEN 

554 FRSTCALL • .FALSE. 

555 C 

556 C COMPUTE THE CONST ANTST FOR KSOIL' S POLYNOMIAL EQUATIONS 

557 C 

558 FOR IDAY-1. NSDAY 

559 A< IDAY ) • -0.85 + 0.733 / ET0< IDAY > 

560 B( IDAY ) ■ 4.97 - 0.661 * ET©( IDAY 1 

561 Ct IDAY > • -8.57 ♦ 1.560 t ET0t IDAY ) 

563 Dt IDAY ) • 4.35 - 0.880 » ET0< IDAY ) 

563 END FOR 

564 END IF 

565 C 

566 PE • ET0< DAY ) 

567 KSTORM • (24 - DUR) / 34.0 

568 C 

569 C IF THE GROUND CELL IS A CANOPY THEN COMPUTE EVATRANSPIRATION 

570 C 30* FROM SURFACE LAYER AND 70* FROM ROOT 20M6 

571 C 

572 IFf TYPENOU .EO. 'CANOPY ' ) THEN 

573 EUAP • PE * KCROP * KSTORM 

574 SFEVAP ■ AMINK SFUC*5 , 0.30*EVAP ) 

575 RTEVAP • AMINK RTUC*95 , 0.70*£VAP ) 

576 SFUC • SFUC - SFEVAP/5 

577 RTUC • RTUC - RTEVAP/95 

578 IF ( SU .GT. 0 ) THEN 

579 SU1NF • AMINK SU , SFEVAP ) 

580 SU • SU - SUINF 

581 SFUC - SFUC + SUXNF/5 

583 END IF 

583 RETURN 

584 END IF 

585 C 

586 C IF THIS CELL IS COVERED BY STANDING UATER THEN COMPUTE 

587 C EVAPORATION OF STANDING UATER 

588 C 

589 KSUEVAP • 1 

590 IF t SU .GT. 0 > THEN 

591 EVAP « PE * KSTORM 

593 SUEUAP • AM1N1 < SU , EVAP ) 
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593 SU • SU - SUEUAP 

594 KSUEUAP • 1 PE - SUEUAP ) / PE 

595 PE • PE - SUEUAP 

596 END IF 

597 C 

598 C COMPUTE EUPORATION FOR BARE SOIL 

599 C 

660 IF ( PE .OT. 6 ) THEM 

601 <1R • ( SFUC - SFUP ) / < SFFC - SFUP ) 

602 IF ( HR .IT. 0 > HR- 0 

603 IF ( HR .CT. 1 ) HR • I 

604 KSOIL • A(DAV) + B(DAV) 2 HR ♦ C(OAV) * 

605 & D(DAY) * HRtHRtHR 

606 IF( KSOIt .LT. 0.05 ) KSOIL • 0.05 

607 IF ( KSOIt .OT. 1.00 ) KSOIL • 1.00 

608 SFEUAP • PE * KSOIt * KSTORN * KSUEUAP 

609 SFEUAP - ArilNlC SFUC*5 , SFEUAP ) 

610 SFUC • SFUC - SFEUAP/S 

611 END IF 

eia c 

613 RETURN 

614 END 

> 


HRJHR ♦ 
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Vt; 


615 C 

616 C 

617 C- 

618 C 

619 C 

626 C 

621 C 

622 C 

623 C 

624 

625 

626 

627 

628 

629 

630 

631 

632 

633 

634 

635 

636 

637 

638 C 

639 C 

640 C 

641 C 

642 C 

643 C 

644 

645 

646 

647 

648 C 

649 C 

650 C 

651 C 

652 

653 

654 

655 

656 

657 

658 C 

659 C 

660 C 

661 C 

662 

663 

664 

665 C 

666 

667 

668 

669 

670 

671 

672 

673 C 

674 C 

675 C 

676 

677 C 

678 C 

679 C 

680 C 

681 
682 

683 

684 

> 



THIS SOUBROUTINE GETS ALL THE AVAILABLE AND NECESSARY 
INFORMATION ABOUT THE REOUSTED DATA BASE CELL. 


SUBROUTINE CELLDATA! ROUSTROU, ROUSTCOL, ELEU, SOIL. CATO, 

& SLOPE. TYPE. CROP, FIELD ) 

IMPLICIT INTEGER ( A - 2 > 

PARAMETER < NC0LM245 , NFIELD*10 ) 

DIMENSION CROPCODE! 10 ), FC ( 14 ) 

REAL SLOPE. 'RESFEET 

CHARACTER*8 TYPE 
LOGICAL FRSTCALL 

COMMON /BUF 1 / SOILS! NCOL ), ElEVS! NCOL ), CATGSt NCOL) 

SPECIAL COMMON BUF1 

COMMON /FILCOD/ FC 

DATA RESFEET/ 328.08 / * 

DATA CURNTROU / 0 /, FRSTCALL/ .TRUE. / 

DATA CROPCODE/ 3. 5. 5, 4, 4, 3, 2, 2, 1, 6/ 

t-ALFALFA. 2-SOYBEAN, 3-UHEAT fc OATS. 4-CORN, S-SORGHUM, 6-PASTURE 


WHEN CALLED FOR THE FIRST TIME OPEN ALL INPUT & OUTPUT FILES 

IF< FRSTCALL ) THEN 
CALL GETFILES 
FRSTCALL • .FALSE. 

END IF 

IF THE ROW UHICH CONTAINS THE REQUESTED CELL IS NOT READ 
YET, READ THE NEXT RECORD OF ALL THREE DATA BASE NAPS. 

IF ( CURNTROU .LT. RQUSTROU ) THEN 

CURNTROU « CURNTROU + 1 

CALL READREC! FC(2>, SOILS, NCOL ) 

CALL READREC! FC(3>, ELEVS, NCOL ) 

CALL READREC! FC!4 ), CATGS, NCOL ) 

END IF 

THE ROU UHICH CONTAINS THE REQUSTED CELL IS IN THE MEMORY 
EXTRACT NECESSARY INFORMATION. 

ELEV - ELEVS! ROUSTCOL ) 

SOIL • SOILS! ROUSTCOL ) 

CATG ■ CATGS! ROUSTCOL ) 

IF! ROUSTCOL .NE. NCOL ) THEN 

NEXTELEV ■ ELEVS! ROUSTCOL+1 ) 

ELSE 

NEXTELEV • ELEVS! ROUSTCOL-1 ) 

END IF 

SLOPE • ATAN! REAL! ELEV - NEXTELEV ) / RESFEET ) 

SLOPE • ABS! SLOPE ) * 57.2958 

DETRMINE THE SOIL TYPE FROM SOIL MAP CODES 

’ SOIL - SOIL / 30 

DETEMINE THE TYPE OF THE CATEGORY ( NON-AGRICULTURAL , 

BRAE SOIL OR CROP). IF IT IS A CROP TYPE FIND CROP. 

IF! ! CATG .GE. 230) .OR. (CATG .LE. 50) ) THEN 
TYPE - 'NONAGRIC' 

ELSE IF! (CATG . LT . 230) .AND. (CATG .GT. 150) ) THEN 
TYPE - ' BARE50IL ' 
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685 ELSE IF< (CATG .GE. 50 ) .AND. (CATG .IE. *50) ) THEN 

686 TYPE * 'CROP 

68? 1C * CATG / 10 - 5 

688 CROP - CRQPCODEC IC ) 

689 FIELD • I100< CATG , (CATG/I0 * 10 J ) ♦ 1 

690 CATG • CATG /10 * 10 

691 END IF 
698 C 

693 C 

694 RETURN 

695 END 

> 
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696 C 
69? C 

698 C CROPDATA 

699 C 

700 C 

?01 C THIS ROUTINE FINDS HORE INFORMATION ABOUT A CELL UHICH IS 

702 C REGISTERED AS CROP. SUCH AS CROP'S PLANTING DATE. CROP 

703 C CONSTANT K, AND ITS DYNAMIC TYPE BASED ON THE CROP CALANDER. 

704 C 

705 C 

706 SUBROUTINE CROPDATAC CROP, FIELD, DATE. KCROP, TYPENOU ) 

707 C 

708 PARAMETER C NCROP-6 , NFIELD-1© , NSTG-5 ) 

709 IMPLICIT INTEGER ( A - 2 ) 

710 DIMENSION PDATESC NCROP.NFIELD ), STAGEDAY( NCROP.NSTG ) 

711 REAL KEQCONSTC NCROP, NSTG. 8), M, A, KCROP 

713 INTEGER FC< 14 > 

713 CHARACTERS TYPENOU 

714 COMMON /FILCOD/ FC 

715 C 

716 DATA <C PDATES ( IC , IF >, IF- 1, NFIELD ), IC-1, NCROP) 

717 C 1-ALFALFA, 2-SOYBEANS, 3-UHEAT & OATS, 4-CORN, 5-SORGHUM, 6-PASTURE 


718 

& 

/ 63. 

66, 

69, 

73. 

71, 

75. 

78, 

82, 

87, 

93. 

719 

& 

126, 

134, 

142. 

149, 

155, 

158, 

162, 

167, 

173, 

186 , 

720 

& 

54. 

59, 

62. 

65, 

68, 

70, 

74. 

79, 

84, 

93, 

721 

& 

107, 

112, 

117, 

122, 

127, 

130, 

135, 

142, 

149, 

161. 

722 

& 

120, 

139, 

134, 

138, 

143. 

149, 

155, 

161, 

169. 

181. 

733 


73. 

73. 

73. 

?3. 

73, 

73, 

73, 

73. 

73. 

73 / 

734 C 




« 








725 


DATA 

(CSTAGEDAYC 

IC, 

IS >. 

IS- 

1, NSTG), 

IC-1, 

NCROP) 

736 

t> 

/ 20, 

61, 

183, 

217, 

365. 






727 

& 

19, 

49, 

102. 

139, 

365. 






728 

& 

31. 

67. 

105. 

117. 

365, 






729 

& 

31, 

73. 

122, 

166, 

365, 






730 

& 

38. 

68. 

113, 

140. 

365. 






731 

& 

20, 

61. 

183, 

317, 

365/ 







733 C 

733 DATA <<C KEQCONSTC IC.IS.IK), IK-1, 2), IS-l.NSTG), IC-1. NCROP) 


734 

& 

/ 0.0 , 

0.7 


0.004, 

0.627, 

0.0 

0 

0.85 

* 

735 

& 

-0.01 , 

2.7 


0.0 . 

0.5 . 





736 

& 

0.0 , 

0.7 


0.018. 

0.478, 

0.0 

t 

1.05 

0 

737 

& 

-0.016. 

2.7 


0.0 , 

0.45 , 





738 

& 

0.0 . 

0.7 


0.011, 

0.356, 

0.0 

0 

l.l 


739 

5 

-0.071, 

8.54 


0.0 , 

0.25 , 





740 

& 

0.0 , 

0.7 


0.010, 

0.398, 

0.0 

0 

1.1 

0 

741 

t. 

-0.013, 

2.625 


0.0 , 

0.55 , 





742 

& 

0.0 , 

0.7 


0.009, 

0 . 455, 

0.0 

0 

1.05 

0 

743 

& 

-0.020, 

3.25 


0.0 , 

0.5 , 





744 

& 

0.0 , 

0.7 


0.001, 

0.68 , 

0.0 

0 

0.75 

0 

745 

& 

-0.007, 

2.1 


0.0 , 

0,5 / 






746 C 

747 C 

748 C. 

749 PLNTDATE • PDATES < CROP , FIELD ) 

750 CROPCNT • DATE - PLNTDATE 

751 C 

753 C DETERMINE THE STAGE OF THIS CROP SUCH AS EMERGED, HARVESTED, .. 

753 C 

754 FOR STAGE-1, 5 

755 IF ( CROPCNT , LT. STAGEDAYC CROP. STAGE ) ) GOTO 10 

756 END FOR 

757 10 M - KEQCONSTC CROP, STAGE, 1 > 

758 A • KEQCONSTC CROP, STAGE, 2 ) 

759 KCROP -AT M * CROPCNT 

760 C 

761 C TREAT ALL AS MEDIUM ROUGH BARE BEFOR EMERGENCE 

762 C 

763 IF C STAGE .LT. 3 ) THEN 

764 TYPENOU • 'MEDMBARE' 

765 RETURN 

> 
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776 C 

777 C 

778 C 

779 

780 

781 
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783 
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END IF 

IF CROP STAGE IS AFTER EMERGENCE AND BEFOR HARVEST 
THEN TREAT IT AS A CANOPY 

IF ( STAGE .17. 5 ) THEN 
TYPENOU • ' CANOPY ' 

- RETURN 
END IF 

AFTER HARVEST TREAT UHEAT, OATS AND ALFALFA AS SMOOTH BARE 
AND TREAT ALL OTHERS AS MEDIUM ROUGH BARE 

IF ( (CROP ,EO. 1 ) .OR. (CROP .£Q. 3) > THEN 
TYPENOU « 'SiTTHBARE' 

ELSE 

TYPENOU • 'MEDMBARE' 

END IF 

RETURN 

END 
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787 C 




788 C 




789 C 

--- 

SOILDATA 


790 C 




791 C 




79a c 


THIS ROUTINE FINDS THE BULK DENSITV OF THE GIUEN SOIL 


793 C 


BASED ON THE BULK OENSITV DISTRIBUTION UITHIN THAT 


794 C 


SOIL, AND RETURNS OTHER REQUIRED INFORMATION ABOUT THE GIUE 


795 C 


SOIL SUCH AS UIITING POINT UATER CONTENT, FIELD CAPACITY 


796 C 


WATER CONTENT, HYDRAULIC CONDUCTIVITY AT SATURATION, ETC. 


797 C 




798 


SUBROUTINE SOILDATAI SOIL. SFPROS, RTPROS, SFFC, RTFC, 


799 

& 

SFWP, SATHC, SUCTION ) 


800 C 




801 C 




80S 


IMPLICIT INTEGER < A - Z ) 


803 C 




804 


PARAMETER ( NSOIL- 8 , NBD- 15 ) 


805 


REAL SFBULK(NBD), B(NSOIL), FSUCTIONCNSQIL ) , ALFA(NSOIL) 


806 


REAL RTBULK ( NBD ), SSUCTION< NSOIL ), SHYDCOND< NSOIL ) 


807 


REAL SFPROS, SATHC, SUCTION, SATSUCT, SFUP. SFFC, RTFC 


808 


REAL RTPROS, FBD, RTBD 


809 


LOGICAL FRSTCALL 


810 


COMMON /BUF2/ PROBABIL< NSOIL , 100 ), PERCENTS ( NSOIL , NBD ) 

811 


SPECIAL COMMON BUF2 


812 C 




813 


DATA ( (PERCENTS! IS, IBD ) , IBD- 1 , NBD), IS-1, NSOIL) 


814 

4 

/ 0, 0, 0, '0. 0, 7, 7, 22, 22, 14, 14, 14, 0, 0, 

0, 

815 

4 

0, 0, 0, 11, 26, 26. 26, 11, 0, 0. 0. 0, 0. 0, 

0. 

816 

4 

2, 3, 3, 8, IB, 20, 17, 11, 12. 3, 3, 0, 0. 0, 

0. 

817 

4 

0, 0, l, 9, 22. 25, 14. It. 9, 5, 2, 2, 0. 0, 

0, 

818 

4 

0, 0, 0, 18, 18, 28, 18, 18, 0, 0, 0, 0, 0, 0, 

0, 

819 

4 

1, 1, 1, 7. 18, 22, 23, 15, 7. 4, 1, 0, 0, 0, 

0, 

820 

4 

0. 0, 0. 7, 19, 18, 30, 15, 7, 4, 0. 0, 0, 0, 

0» 

821 

4 

0, 0, 0, 0. 0, 6, 0, 0, 0, 15, 31, 23, 8, 15, 

8 / 

822 C 




823 


DATA SFBULK/ 0.818, 0.888, 0.957, 1.03, 1.10, 1.17, 1.24, 

1.30 

824 

4 

1.37 , 1.44 , 1.51 , 1.58. 1.65, 1.72. 1.79 

/ 

825 • 


DATA RTBULK/ 1.10, 1.15, 1.20, 1.25, 1.30, 1.35, 1.40. 1.45. 

826 

4 

1.50, 1.55, 1.60, 1.65, 1.70, 1.75, 1.30 / 


827 C 




828 


DATA B 


829 

4 

/ 9.77, 6. 66, 7.21, 6.31, 6.22, 6.81, 6.66, 4.26 / 


830 


DATA FSUCTION 


331 

4 

/ 2. .49, 13.91. 12.90, 13.18, 16.03, 13.18, 13.91, 3.22 / 


832 


DATA SSUCTION 


833 

4 

/ 31.20. 34.67, 30.76, 33.08, 39.25, 33.08, 34.67, 10.14 / 


834 


DATA SHYDCOND 


835 

4 

/ 0.814, 2.33, 2.02,. 2.31, 2.39, 2.31, 2.33, 58.36 / 


836 


DATA FRSTCALL/ .TRUE. / 


837 C 




838 C 




839 C 


SET 1JP THE PROBABILITY ARRAY FOR BULKDENSITY DETERMINATION 


840 C 


OF EACH SOIL TYPE. (ONLY AT FIRST CALL) 


841 C 




842 


IF< FRSTCALL ) THEN 


843 


FRSTCALL • .FALSE. 


844 


FOR IS-1, NSOIL 


845 


ALFA( IS ) • 1 / B( IS ) 


846 


END FOR. 


847 


FOR IS-1, NSOIL 


848 


START • 1 


849 


STOP ■ 0 


850 


FOR BDCODE-l , NBD 


851 


PERC ■ PERCENTSt IS . BDCODE ) 


852 


IF ( PERC .NE. 0 ) THEN 


853 


STOP ■ STOP + PERC 


854 


FOR I -START, STOP 


855 


PROBABILt IS, I ) • BDCODE 


856 


END FOR 



> 
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85? START • START ♦ PERC 

858 END IF 

859 END FOR 

860 IF< STOP .HE. 100 ) URITEM 1, ' < IX# ‘STOP- *,14, 16 ) ' JSTOP, IS 

861 END FOR 

862 END IF 

863 C 

864 C GET THE BULK DEN5ITV ACCORDING TO THE PROBABILTV FOR THIS SOIL 

865 C FOR BOTH SURFACE LAVER <0-6 CPU AND ROOT LAVER <5-95 CPU 

866 C 

86? RANDOM -IRAN< 1 , 100 > 

868 BDCODE ■ PROBABILf SOIL, RANDOPI ) 

869 SF3D • SF0ULK ( BDCODE ) 

870 RTBD» RTBUUCC BDCODE > 

8?1 C 

872 C CALCULATE THE UATER CONTENT OF EACH SOIL TVPE AT UILTINO POINT 

873 C AND FIELD CAPACITV. SPECIFIC BULK DENSITV IS 2.65 FOR ALL SOILS 

874 C 

875 SFPROS ■ 1 - SFBD / 2.65 

876 RTPROS ■ 1 - RTBD / 2.65 

87? SFFC - SFPROS * < SSUCTION< SOIL ) / 333.0 ) tt ALFA< SOIL ) 

878 RTFC • RTPROS t < 5SUCTI0NC SOIL ) / 333.0 ) tt ALFA< SOIL ) 

879 SFUP ■ SFPROS * < 5SUCTI0N< SOIL ) / 15000 ) Xt ALFA< SOIL ) 

880 SUCTION • FSUCTIOfU SOIL > 

881 SATHC • SHVDCONDC SOIL ) 

882 RETURN 

883 END 

> 
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884 C 

885 C 

886 C 

88? C 

888 C 

889 C 

890 C 

891 C 

892 

893 C 

894 

895 

896 

897 

898 C 

899 C 

900 

901 

902 

903 

904 

905 

906 

907 

908 
989 

910 

911 

912 

913 

914 

915 C 

916 

917 

918 101 

919 102 

920 103 

921 104 

922 

923 


GETFILES 

THIS ROUTINE OPENS ALL THE INPUT AND OUTPUT FILES AND 
ASSIGNS AN AVAILABLE UNIT (FILECODE) ON UHICH THE FILE 
(JILL BE OPENED. 

SUBROUTINE GETFILES 

INTEGERS FILENAME! 17 1 
INTEGER FC! 14 >, ERRCODE 
LOGICAL ERR 
common /FILCGD/ FC 


URIT£( 11, 104 ) 

READ< 12, 102 > FILENAME 

CALL OPN! FC(1), FILENAME, 'OLD ' , 'FOR ' , ERRCODE, ERR > 
IF ( ERR 1 THEN 

uRiTEUi.ien filename, errcode 

STOP 
END IF 

READCFC! 1 ) , * ) NUMFILES 
FOR I »2,NUMFILESM 
R£AD!FC< 1 ;, 102 ) FILENAME 

CALL OPN(FC! I ) ,FILENAMS, 'OLD' , 'UNF' .ERRCODE, ERR ) 

IF t ERR > THEN 

WRITE! 11. 101) FILENAME, ERRCODE 

STOP 

END IF 

WRITE! 11, 103 1 FILENAME, FClI) 

END FOR 


RETURN 

FORMAT! IX, 'ERROR S*SS UHILE 
FORMAT! 17A1 ) 

FORMAT (1 X. 17A1, 'UAS ASSIGNED 
FORMAT! IX, 'ENTER NAME OF THE 
fc ' AND OUTPUT FILE NAMES') 


QPENINIG ' 17A1, 'ERRCODE* ' 13 ) 
TO'24) 

FILE UHICH CONTAINS INPUT', 


END 
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924 C 

925 C 

926 C- 
92? C 

928 C 

929 C 

930 C 

931 C 

932 

933 C 

934 

935 

936 

937 

938 

939 C 

940 C 

941 C 

942 

943 C 

944 

945 

946 

947 

948 

949 

950 

951 C 

952 C 

953 

954 C 

955 C 

956 

957 

958 


959 

960 

961 

962 

963 

964 

965 

966 
96? 

968 

969 

970 

971 C 

972 101 

973 102 

974 

> 
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THIS ROUTINE HAS TUO ENTRIES USED IN READING AND URITING 
FROn AND TO THE I/O FILES 

SUBROUTINE I0CALLS 

IMPLICIT INTEGER t A - Z ) 

PARAMETER ( NURD -1245 , NPAUSE-3 ) 

DIMENSION RECORD< NURD >, MATRIX! NPAUSE, NURD ) 

INTEGER PC ( 14 ) 

COMMON /FILCOD/ FC 


ENTRY READREC ( FILECODE, RECORD, NUORDS ) 


BUFFER IN < FILECODE. RECORD, B, NUORDS, 10 ) 
CALL STATUS! FILECODE ) 

IF! 10 .NE. 2 > THEN 

URITE (11,101 ) FILECODE, 10 

STOP 

END IF 

RETURN 


ENTRY URITDATA ( IDENT , MATRIX, NCOL > ' 


IF! IDENT .EO. 1 ) 0TFILE*FC(6) 

IF! IDENT .EQ. 2 ) OTFILE-FC! 10 ) 

FOR IP-1, NPAUSE 
FOR COL- l, NCOL 

RECORD! COL ) • MATRIX! IP , COL ) 

END FOR 

BUFFER OUT! OTFILE , RECORD, B, NCOL, 10 ) 

CALL STATUS! OTFILE ) 

IF! 10 .NE. 2 ) THEN 
URITE! 11, 102 ) OTFILE, 10 
STOP 
END IF 

OTFILE • OTFILE + 1 

END FOR 

RETURN 

FORMAT! IX, 'ERROR tXXX UHILE READING FROM UNIT'14,' STATUS', 14) 
FORMAT! IX, 'ERROR XXtX UHILE URITING TO UNIT '14,' STATUS-'I4) 

END 
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975 C 

976 C 

977 C UPDTHIST 

978 C 

979 C 

980 C THIS ROUTINE 'UPDATES THE HISTOGRAMS 

981 C 

98S SUBROUTINE UPDTHISTC UAL, HIST, INDEX, MIN, MAX ) 

983 IMPLICIT INTEGER ( A - 2 ) 

984 PARAMETER ( NPAUSE-3 ) 

985 DIMENSION HIST< NPAUSE, MAX ) 

986 C 

987 C 

988 C UPDATE MIN, MAX UALUES OF DATA 

989 C 

990 HIST < INDEX . 1 ) - MIN0( HISTC INDEX . 1 > . UAL ) 

991 HIST < INDEX , MAX-MIN+S) « MAX0< HIST< INDEX . MAX-NIN*5>, UAL ) 

992 C 

993 C UPDATE THE TOTAL NUMBER OF UALUES COUNTED 

994 C 

995 HIST< INDEX , HAX-MIN+6 ) • HIST( INDEX, MAX-MIN+6 )♦! 

996 C 

997 C UPDATE THE FREQUENCY COUNT FOR THIS UALUE 

998 C 

999 UAL - MAX0< MIN0C MAX*1 , UAL), MIN-l ) 

1000 HIST( INDEX, UAL-MIN+3 ) - HIST( INDEX. UAL-MIN+3 ) ♦ 1 

1001 RETURN 

1008 END 

> 
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1003 C 

1004 C 

1005 C- 

1006 C 
100? C 

1008 C 

1009 C 

1010 C 

1011 C 
1018 

1013 C 

1014 

1015 

1016 

1017 C 

1018 C 

1019 
10S0 

1021 C 

1022 C 

1023 C 

1024 C 

1025 

1026 C 
102? C 

1028 C 

1029 C 

1030 

1031 

1032 

1033 

1034 

1035 

1036 

1037 C 

1038 C 

1039 C 

1040 C 

1041 

1042 

1043 

1044 C 

1045 C 

1046 C 

1047 C 


1048 

1049 

1050 

C 


1051 

100 


1052 


& 

1053 

101 


1054 

102 


1055 

103 


1056 


8. 

1057 
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OTPTHIST 


THIS ROUTINE TAKES A GIVEN HISTOGRAM ARRAY CONTAINING THE 
FREQUENCIES AND WRITES OUT THE PERCENTS AND CUMULATIVE 
PERCENTS FOR EACH DATA VALUE. 

SUBROUTINE OTPTHIST! HIST, INDEX, NXN, MAX > 

IMPLICIT INTEGER ( A - Z ) 

DIMENSION HIST ( 3, MAX > 

REAL TOTAL, SUM 


TOTAL • HIST ( INDEX , MAX-MIN+6 > 

SUM • HIST! INDEX , 2 ) 

URITE OUT INFORMATION ON DATA POINTS ENCOUNTERED UHICH 
WERE LESS THAN THE INDICATED MINIMUM VALUE 

URITEU3.100) MIN, HIST! INDEX, 2 ), SUM/TOTAL. SUM/TOTAL 

URITE OUT PERCENTS AND CUM PERCENTS FOR ALL UAlUES FROM 
INDICATED M1NMUM VALUE THROUGH MAXIMUM VALUE 

FOR PTR- 3, MAX-MIN»3 
COUNT • HIST ( INDEX. PTR > 

SUM * SUM ♦ COUNT 

IF 1 COUNT .HE. 0 ) THEN 

URITE( 13, 101 ) MIN+PTR-3, COUNT, COUNT'TOTAL, SUIVTOTAL 
END IF 
END FOR 

URITE OUT INFORMATION ON DATA POINTS ENCOUNTERED UHICH 
UERE LARGER THAN THE INDICATED MINIMUM VALUE 

COUNT • HISTt INDEX, MAX-MIN+4 ) 

SUM • SUM + COUNT 

URITE! 13, 102 > MAX, COUNT, COUNT/TOTAL, SUM/’TOTAL 

URITE OUT TOTAL NUMBER OF VALUES, MIN VALUE AND MAX VALUE 
THAT UAS ENCOUNTERED 

URITE! 13, 103 ) TOTAL, HIST< INDEX , 1 >, HIST< INDEX , MAX-MIN+S) 
RETURN 

FORMAT U9X, 'RANGE' , 14X' COUNT' ,?X, 'PERCENT' 10X'CUM PERCENTS' . 
/✓5X, ' < ' 17, 5X, '*•' ,5X, I8,5XE14 . 7, SXE14.7 ) 

FORMAT (7X, I6,SX,'-'SX, I8.5XE14 . 7, 5X.E14. 7 ) 

F0RMAT(5X, ' > ' , 17, SX, '-',SX, I8.SXE14.7.5XE14.7) 

FORMAT!2XF10.1 , ' TOTAL VALUES COUNTED'SX'MIN AND MAX VALUES' 

' ENCOUNTERED • '19,2X19) 

END 


NASA-JSC 


127 


si 


V 


